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ABSTRACT 

The  non-linear  elasto-plastic  response  of  a  clamped 
flat  plate  to  a  typical  air  shock  wave  were  investigated. 
The  nonlinear  effects  to  the  plate  responses  due  to  the 
material  and  geometric  nonlinearity  were  studied. 

In  this  study ,  (1)  the  necessity  of  the  modification  of 
old  armored  vehicles  was  reviewed  and  (2)  the  NASTRAN  code 
was  employed  in  this  investigation,  (3)  the  theoretical 
background  of  the  nonlinear  transient  analysis  was  described, 

(4)  a  step  by  step  procedure  of  analyzing  the  dynamic  load 
problem  by  shock  wave  using  NASTRAN  code  was  developed,  and 

(5)  sensitivity  analyses  were  performed  and  also  difficulties 
associated  with  the  nonlinear  analysis  were  described. 
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I.   INTRODUCTION 

Conventional  tanks  (old  tanks)  which  have  steel  armor, 
can  be  completely  defeated  by  one  round  of  modern  projectiles 
such  as  anti-tank  missiles.   That  has  been  proved  in  the  1973 
Arab-Israeli  war.   In  that  war,  the  40%  of  total  Israeli 
casualties  were  caused  by  armored  vehicle  destructions.   Thus, 
tank  armor  improvement  became  the  first  priority  in  the  Israeli 
tank  corps . 

If  war  occurs  in  any  country,  the  same  problem  will  be 
faced.   That  is  the  reason  why  the  modification  of  the  conven- 
tional tank  armor  is  needed  against  modern  anti-tank  weapons . 
To  improve  the  armor  protection,  the  material  properties, 
projectile  characteristics,  and  explosion  phenomena  must  be 
understood  well.   This  is  a  complicated  and  difficult  problem 
and  is  originally  a  nonlinear  problem. 

First  the  shock  wave  from  an  explosion  imposes  a  dynamic 
load  on  any  object.   That  dynamic  load  is  characterized  by 
rapidly  reached  peak  overpressure  which  then  decreases  as  the 
shock  wave  decays.   The  net  effects  of  the  load  depend  both 
on  the  nature  of  the  shock  wave  and  the  geometry  and  construc- 
tion of  the  object.   For  example,  if  the  explosion  occurs  on 
thin  skirting  plate,  perforation  takes  place  by  impact  load 
and  their  impulse  loading  for  main  armor  would  be  reduced 
significantly.   The  typical  air  shock  wave  was  characterized 
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as  follows;  the  shock  wave  is  supersonic  and  generates  the 
abrupt  increasing  overpressure. 

The  non-linear  elasto-plastic  responses  of  a  clamped  flat 
plate  to  a  typical  air  shock  wave  were  investigated.   The 
nonlinear  effects  to  the  plate  responses  due  to  the  material 
and  geometric  nonlinearity  were  studied. 

In  this  study,  (1)  the  necessity  of  the  modification  of 
old  armored  vehicles  was  reviewed,  (2)  the  NASTRAN  code  was 
employed,  (3)  the  theoretical  background  of  the  nonlinear 
transient  analysis  was  described,  (4)  a  step  by  step  proce- 
dure of  analyzing  the  dynamic  load  problem  by  shock  wave 
using  NASTRAN  code  was  developed,  and  (5)  sensitivity  analy- 
ses were  performed  and  also  difficulties  associated  with  the 
nonlinear  analysis  were  described. 


II .   LITERATURE  REVIEW 

A.   MODERN  PROJECTILES 

An  important  version  of  the  AP  (Armor-Piercing)  round  is 
that  known  as  APDS  (Armour  Piercing  Discarding  Sabot) .   In 
this  round  the  armour-piercing  element  is  a  hard  core  of 
calibre  significantly  smaller  than  that  of  the  complete 
projectile  and  shaped  for  maximum  anti-armour  effectiveness. 
This  is  surrounded  by  a  much  lighter  casing,  or  sabot,  which 
breaks  up  as  the  projectile  leaves  the  gun  barrel.   The 
point  of  this  arrangement  is  that  the  core  can  have  a  shape 
which  differs  considerably  from  that  of  the  complete  projec- 
tile and,  not  having  been  scored  by  the  rifling  in  the  barrel, 
is  aerodynamically  clean.   Moreover,  since  it  is  very  much 
heavier  than  the  sabot,  the  core  acquire  a  greater  energy 
than  could  be  imparted  to  it  if  it  were  fired,  without  sabot, 
by  a  gun  of  appropriate  calibre  and  the  same  chamber  pressure 
as  that  of  the  larger  gun.   A  controversial  refinement 
involves  the  use  of  depleted  uranium  in  the  core.   This 
metal  has  an  exothermic  reaction  with  steel  which  aids 
penetration;  it  is  controversial  because  of  its  connection 
with  the  other  isotopes  of  uranium  and  their  radioactive 
properties. 

High  explosive  (HE)  rounds  with  contact  fuzes  can  also 
be  used  against  tanks  and  may,  if  they  are  large  enough,  be 
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effective  against  exposed  portions  of  the  tank  mechanism 
(such  as  the  tracks)  but  will  not  in  general  be  effective 
against  heavy  armour.   Two  modern  refinements  of  the  HE 
round  may  be  effectively  used  in  some  circumstances:   they 
are  the  HEAP  (High-Explosive  Armour  Piercing)  and  HESH 
(High-Explosive  Squash-Head)  rounds. 

The  HEAP  round  has  a  projectile  of  relatively  high 
velocity,  a  hardened  case  and  a  high-explosive  charge  which 
is  contact  fuzed  with  a  very  short  delay.   The  effect  of  the 
arrangement  is  that  the  projectile  either  embeds  itself  in 
the  armour  or  passes  through  it  and  then  explodes  with  an 
effect  which  is  more  damaging  than  that  of  the  simple  HE 
round . 

HESH  rounds  can  be  used  with  lower-velocity  weapons,  and 
are  so  designed  that  the  nose  of  the  shell  crumples  on  impact 
to  make  a  large  contact  area  before  the  explosive  is  detonated 
Such  rounds  will  not,  in  general,  pierce  armour,  but  the 
explosion  sets  up  shock  waves  in  the  armour  which  in  them- 
selves can  incapacitate  the  crew  and  may  also  cause  pieces  of 
metal  to  break  away  from  the  surface  remote  from  the  explosion 
Such  fragments  fly  off  at  high  velocity  and  can  do  great 
damage  inside  a  tank  turret. 

Most  widely  used  in  infantry  weapons,  however,  is  the 
projectile  (which  may  be  a  shell,  a  rocket,  a  grenade  or  a 
guided  missile)  with  a  shaped-charge  warhead  frequently 
described  as  a  HEAT  (High  Explosive  Anti  Tank)  round.   This 
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type  of  projectile  is  designed  to  be  effective  with  only  a 
moderate  impact  velocity  and  depends  for  its  effect  on  the 
design  of  the  warhead  and  is  formed  round  a  hollow  cone  with 
its  apex  to  the  rear.   The  base  of  the  cone  is  of  approxi- 
mately the  same  diameter  as  the  warhead  at  that  point  and  in 
front  of  it  is  an  air  space  enclosed  by  the  rounded  or 
pointed  head  of  the  projectile  on  which  the  fuze  igniter  is 
mounted.   The  operation  of  the  fuze  is  critical  to  the  success 
of  the  charge.   It  has  to  have  a  set  distance  from  the  surface 
of  the  armour  when  it  is  detonated  and  for  this  reason  the 
igniter  is  carried  well  in  front,  often  giving  rise  to  the 
distinctive  pointed  nose  of  a  hollow  charge  round. 
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Figure  2.1.   The  Effects   of  Current  Anti-Tank  Projectiles 

on  the  Steel. 


The  igniter  is  connected  to  the  detonator  which  is  at  the 
back  of  the  charge,  this  starts  a  detonation  wave  running 
through  the  conical  explosive  and  the  result  is  that  the 
force  of  the  explosion  is  concentrated  into  a  thin  jet 
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moving  at  a  speed  of  oyer  20,000  m/s  and  at  a  great  tempe- 
rature.  The  jet  takes  a  finite  distance  to  form,  hence  the 
importance  of  having  the  igniter  in  front  and  the  distance 
required  to  form  a  jet  is  known  as  the  "stand  off".   For  most 
modern  warheads  the  stand  off  distance  needs  to  be  between  2 
and  3  diameters  of  the  warhead  and  to  provide  the  full  dis- 
tance sometimes  interferes  with  the  ballistics  of  the 
projectile . 

The  effects  of  a  hollow  charge  on  the  interior  of  a 
vehicle  can  be  variable.   A  very  small  jet  may  well  not 
cause  much  damage  apart  from  making  a  hole  in  the  armour  and 
destroying  whatever  is  in  its  direct  path,  but  it  is  possible 
to  "turn"  warheads  so  that  the  penetrative  effect  is  lessened 
but  there  is  more  energy  left  after  passing  through  the 
armour  and  this  energy  can  be  spread  out  into  more  of  a  fan- 
shaped  exit  blast  with  appropriately  enhanced  damage.   All 
hollow  charge  jets  carry  a  certain  amount  of  molten  metal 
into  the  interior  and  some  research  has  gone  into  finding 
ways  of  carrying  more  through  the  hole. 

Two  examples  of  the  Soviet  anti-tank  guided  missiles  are 
described  below.   [Ref.  1] 

1.   Sagger  Anti-Tank  Guided  Missile 

Sagger  is  the  NATO  code-name  for  the  Soviet  anti- 
tank missile  which  is  also  known  by  the  US  alphanumeric 
designation  AT-3,  which  is  now  extensively  deployed  in 
Europe  and  elsewhere.   Currently  it  is  known  to  be  used  in 

at  least  three  ways: 
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a)  As  a  man-portable  missile 

b)  Mounted  on  a  bracket 

c)  Mounted  on  a  vehicle 

The  characteristics  of  sagger  missile  is  as  follow. 

a)  Type:   Wire-guided  MCLOS  (Manual  Command  to  Line 
of  Sight) 

b)  Weight  at  launch:   11.3  kg 

c)  Diameter:   120  mm 

d)  Warhead:   HEAT 

e)  Range:   500-3000  m 

f)  Penetration:   400  mm 

2 .   AT-4  Anti-Tank  Guided  Missile 

There  have  been  reports  for  some  years  of  a  Soviet 
SACLOS  (Semi -Automatic  Command  to  Line  of  Sight)  system,  but 
it  was  not  until  1980  that  photographs  and  a  little  informa- 
tion appeared  in  some  Warsaw  Pact  military  journals. 

Characteristics : 

a)  Type:   Wire  guided  SACLOS  system 

b)  Weight  at  launch:   40.4  kg 

c)  Diameter:   120  mm 

d)  Warhead:   HEAT 

e)  Range:   2000  -  2500  m 

f)  Penetration:   500-600  mm 
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Figure  2.2.   Sectioned  View  of  Sagger  Missile. 

Which  of  these  various  types  of  projectile  will  be  most 
effective  against  armour  which  depends  on  the  design  and 
thickness  of  the  armour.   Against  solid  armour  the  HEAT 
round  is  particularly  effective;  but  in  certain  circumstances 
a  sandwich  construction  is  likely  to  be  more  vulnerable  to  an 
APDS  round.   There  is  therefore  a  good  deal  to  be  said  for 
having  a  mix  of  projectiles  available  for  satisfactory 
defense  against  tanks. 

B.   CONFIGURATION  OF  ARMORED  VEHICLES 
1.   Armor 

We  can  define  the  military  armor  as  a  material 
intended  to  protect  men  and  installations  of  equipment. 
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Aluminum,  titanium  and  magnesium  alloy,  plastic  laminates, 
nylon  fabrics  and  masticstone  composites  have  all  been  used, 
but  the  most  common  form  of  armor  is  heat-treated  alloy  steel. 
This  comes  mainly  in  two  forms:   homogeneous  and  face-hardened 

The  first  has  constant  ballistic  and  mechanical  quali- 
ties throughout,  which  in  modern  tank  armor  means  high  tensile 
strength  (of  the  order  of  150,000  lb  sq  in)  and  medium  hard- 
ness (i.e.  Brinnell  N  300);  face-hardened  armor  consists  of 
a  very  hard  first  layer  which  resists  penetration  but  lacks 
toughness  and  a  more  elastic  core  to  support  the  face.   Naval 
deck  armor  and  modern  tank  armor  are  homogeneous,  while  face- 
hardened  armor  is  still  used  for  some  light  armored  vehicles 
and  side  protection  for  naval  vessels. 

Steel  armour  usually  means  rolled  or  wrought  plates, 
but  castings  are  used  for  tank  turrets  and  even  for  entire 
tank  hulls,  as  in  the  M.48  and  M.60  tanks  (Steel  armor  is 
usually  a  nickel-chrome-molybdenum  alloy,  but  other  alloys 
including  manganese-silicon-chromium-nickel  steel  have  been 
used,  especially  in  Soviet  tank  armor.) . 

Typically,  a  well-designed  hollow-charge  warhead  can 
be  expected  to  penetrate  armor  up  to  five  times  the  diameter 
of  the  warhead.   The  best  defense  against  such  attack  is  not 
to  make  thicker  armor,  but  to  dissipate  the  jet  before  it 
reaches  the  surface.   To  do  this  a  light  screen  is  carried 
some  distance  in  front  of  the  main  armor,  of  sufficient 
thickness  to  set  off  the  fuze  of  any  hollow-charge  projectile. 
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Quite  naturally  such  screens  do  not  last  long,  but  they  are 
easily  replaced  and  provide  almost  complete  protection  for 
the  armor  behind.   A  common  use  of  such  a  technique  is  the 
thin   skirting  plates  carried  over  the  tracks  of  all  modern 
tanks.   Before  leaving  this  subject  it  is  appropriate  to 
mention  a  recent  British  armor  development,  the  'Chobham' 
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armor.   This  is  claimed  to  offer  greatly  enhanced  resistance 
to  all  forms  of  projectile  as  is  indicated  by  the  accompanying 
illustrations.   Details  of  the  construction  of  this  armor  have 
not  been  released,  but  it  seems  probable  that  it  is  a  sand- 
wich made  up  of  two  sheets  of  armor  plate  separated  by  a 
highly  dispersive,  possibly  granular,  filling.   The  disad- 
vantage of  this  armor  lies  in  its  bulk.   To  gain  the  protec- 
tion, it  has  to  be  thick  and  fairly  heavy;  tanks  are  already 
large  enough  and  the  addition  of  armor  plating,  which  is  200 
mm  or  more  in  thickness  makes  a  startling  difference  to  the 
outline  and  bulk  of  the  entire  vehicle. 
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2 .   Current  Modification  of  Old  Tanks 

Determining  which  one  should  we  choose,  a  new  one 
or  modification  of  an  old  one  usually  depends  on  cost  and 
political  situations. 

The  U.S.  completed  conversion  of  M48  tanks  to  M48A5 
standards  in  19  80.   This  conversion  basically  included  adding 
diesel  engines  to  M48Als  and  adding  105-mm  guns  to  the  tur- 
rets of  M48Als  and  M4  8A3s.   Numerous  minor  changes  were  also 
made.   But  little  improvement  in  armor  protection  has  been 
made.   If  we  want  to  take  more  modifications  for  M4  8A3s  and 
M4  8A5s  tanks,  a  spaced  applique  armor  should  be  added  and 
include  a  1,20  0  hp  engine,  improved  transmission,  better 
fire  control,  and  a  hydropneumatic  suspension  system. 

The  following  example  shows  modern  modifications. 

The  High  Performance  M60  MBT  has  been  developed  as  a 
private  venture  by  the  General  Products  Division  of  Teledyne 
Continental  Motors.   It  is  essentially  an  M60  series  MBT 
fitted  with  a  new  powerpack,  new  suspension  and,  as  an  option, 
additional  armor.   These  modifications  can  be  carried  out  on 
existing  M60  series.   Teledyne  believes  that  to  fit  an  M60A1 
with  the  new  powerpack  would  cost  under  $200,000,  while  the 
added  armor  would  cost  just  $40,000. 
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SPACED  APPLIQUE  ARMOR 
EQUALS  40  64  CM   RHA 


SPALL  ■> 


<W!M 


Figure  2.4.   Applique  Armor  in  M60  MBT . 

3 .   Israel  Add-On  Armor 

Israeli  experience  in  the  196  7  campaign  proved  that 
mobility  was  not  substitute  for  armor  protection  and  they 
therefore  decided  at  an  early  stage  that  the  main  emphasis 
would  be  placed  on  armor,  with  firepower  and  mobility  second 
and  third  priorities.   After  that  they  spent  a  lot  of  money 
to  figure  out  this  protection  problem.   Total  cost  of 
research,  development,  trials  and  the  building  of  the  proto- 
type vehicles  was  about  $65  million. 

In  May  1977  Israel  finally  announced  that  it  had 
developed  a  new  MBT  called  the  Merkava  which  was  made  of 


19 


cast  and  welded  armor  with  a  well-shaped  glacis  plate. 
Behind  the  first  layer  of  cast  armor  is  a  space  filled 
with  diesel  fuel  and  then  another  layer  of  armor.   This 
spaced  armor  gives  the  tank  protection  from  HEAT  projec- 
tiles and  ATGWs  (anti-tank  guided  weapons) . 

They  also  modified  their  old  tanks  such  as  M60s, 
Centurions,  etc. 

The  add-on  and  readily  detachable  protective  arrays 
seen  on  Israeli  M60s  and  Centurions  contain  ceramic  tiles, 
would  also  mark  a  first  in  terms  of  operational  deployment. 
The  possibility  cannot  be  excluded  that  armor  of  both  kinds 
is  involved. 

While  no  answers  are  to  be  found  in  external  appear- 
ance, the  configurations,  as  seen  in  the  accompanying  Defence 
Attache  photographs,  show  obvious  differences  in  the  thick- 
ness and  shape  of  the  packs  provided  for  M60  and  Centurion 
respectively.   The  Centurion's  additional  armor  is  mainly  in 
the  form  of  quite  shallow  panels  fitted  to  the  brow  of  the 
turret  and  thicker  pieces  to  the  glacis  and  the  mantlet 
either  side  of  the  main  gun;  a  hole  is  included  in  the  left- 
hand  piece  to  permit  the  firing  of  the  coaxial  machine  gun. 
Wedge-shaped  packs  are  mounted  on  top  of  the  storage  boxes 
fitted  forward  above  the  track  guards. 
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Figure  2.5.   Add  on  Armor  Blocks  in  M60. 
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Figure  2.6.   Add  on  Armor  Blocks  in  Centurion. 

The  M60,  partly  because  it  lacks  whatever  protective 
values  come  with  the  Centurion's  storage  bins,  has  been 
provided  with  a  considerable  amount  of  add-on  armor  around 
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the  turret.   Here  the  packs  are  thick  and  vertically  deep, 
with  smaller  and  thinner  sections  patched  in  --  with  no 
great  concern  for  aesthetics  --  to  the  upper  surfaces  which, 
because  of  the  turret's  curvature  and  varying  angles  to  the 
perpendicular,  not  to  mention  obstructions  such  as  lifting 
eyes,  cannot  be  covered  by  the  main  boxes.   As  with  the 
Centurion,  panels  are  fitted  to  the  front  of  the  turret  roof 
and  to  the  glacis;  in  the  latter  case,  because  the  M60's 
basic  glacis  is  curved,  it  has  been  necessary  to  fit  a 
straightening  frame  to  which  the  packs  can  be  attached.   The 
M60's  installation  is  completed  by  the  fitting  of  wedge- 
shaped  or  angled  sections  either  side  of  the  driver's  hatch 
to  help  protect  the  vulnerable  area  where  turret  meets  hull, 
with  further  angled  sections  each  side  of,  and  thinner  plates 
above,  the  mantlet. 

Although  it  can  be  stated  with  a  fair  degree  of 
certainty  that  the  add-on  arrays  are  intended  to  help  defeat 
shaped-charge  high-explosive  anti-tank  (HEAT)  attack,  any 
discussion  of  what  the  packs  and  panels  consist  of  must,  in 
the  absence  of  any  word  of  guidance  from  the  Israelis  -  who 
on  the  contrary  are  happily  watching,  if  not  promoting,  the 
circulation  of  a  number  of  often  fanciful  and  frequently 
mutually  exclusive  'solutions'  to  the  puzzle-fall  short  of 
a  definitive  answer. 

While  one  source  suggested  that  the  boxes  were  empty 
(and  the  large  ones  on  the  M60  do  sound  quite  hollow  when 
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thumped)  the  probability  is  that  there  is  more  to  it  than 
that.   In  view  of  solutions  that  have  been  adopted  else- 
where, though  usually  in  the  context  of  new-build  tanks,  it 
seems  quite  likely  that  the  Israelis  have  developed  their  own 
ceramic  tiles  and  that  the  boxes  are  partially  filled  with 
them.   The  alternative  of  a  form  of  active  armor  is  suggested 
more  by  the  vehemence  with  which  several  sources  in  Israel 
stated  this  to  be  the  case  than  by  any  objective  evidence. 

A  number  of  different  types  and  configurations  of 
active  armor  are  possible.   An  example  is  a  sheet  of  explosive- 
laid  of  top  of  a  steel  plate,  with  a  covering  for  environmental 
protection  on  top  -  which  detonates  when  hit  by  the  incipient 
jet  of  a  shaped-charge  warhead  and  in  the  process  disrupts  the 
jet's  further  formation.   If  the  Israelis  are  indeed  deploying 
a  form  of  active  armor,  such  a  fact  would  carry  the  implica- 
tion that  they  have  sufficiently  overcome  the  formidable 
difficulties  that  seem  to  have  prevented  larger  and  better- 
endowed  countries  from  doing  likewise;  the  United  States  and 
Soviet  Union  are  among  the  nations  that  are  known  to  have  been 
engaged  in  related  research  and  experimentation  for  many  years. 

C.   ELEMENTARY  DETONATION  THEORY 

1.   The  following  assumption  will  be  taken  for  the 
simplest  analysis  of  a  detonation:  [Ref .  2  ] 
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a)  The  flow  (motion  of  the  shock  front)  is  one 
dimensional 

b)  The  plane  detonation  front  is  a  jump  discontinuity 
at  which  chemical  reaction  is  assumed  to  be 
instantaneous . 

c)  The  discontinuity  is  time-independent,  so  that  the 
product  composition  (at  least  directly  behind  the 
shock,  front)  is  likewise  time-independent. 

u   represents  the  particle  velocity  of  blast  wind,  which  is 
that  of  the  products  behind  the  shock  front.   Subscript   x 
represents  the  initial  values  of  the  properties,   y   those 
after  shock.   D   is  the  detonation  speed  and   p   is  the  shock 
wave  overpressure. 

From  mass  conservation  across  the  front  we  can  write 


p   D  =  p   (D-u  ) 
x      y     p 


(1) 


Vpy 


Uf 


Vpx 


Figure  2.7.   Explosive  Shock  Front. 
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and  from  conservation  of  momentum, 

p   -  p   =  p   u  D  (2) 

■y   rx    x  p 

Eliminating  up  from  the  two  equations  we  obtain  the  equation 

for  the  Rayleigh  line: 

P  2D2  -  (p   -  p  )/(v   -  v  )  =  0  (3) 

x       ^y    ^x'  /   x     y  v  ' 

(Observe  that  in  part  of  the  previous  equation  we  replace   1/p 

by  the  specific  volume,  v) .   Since  for  a  detonation,  p   >  p  , 

y    x 

it  follows  that   v   <  v   (for  real  values  of  the  detonation 

velocity) ;  and,  indeed,  it  is  observed  as  a  rule  that  there 

is  about  a  one-third  increase  in  density  on  detonation.   If 

this  condition  is  substituted  in  eq  (1),  it  follows  that 

0  <  u   <  d,   so  that  the  products  travel  toward  the  detonation 
p  * 

front,  though  trailing  behind. 

If  the  Rayleigh   equation  is  solved  for   p    we  obtain 
p   =  (p   +  p  D2)  -  Px2D2  v  (4) 

wherein   p    is  linear  with   v  ,  the  slope  of  the  line  being 

2  2 

-P   D  .   The  conservation  conditions  thus  require  that  the 

initial  and  final  states  of  a  detonation  must  lie  at  the 
termini  of  the  Rayleigh  line. 

Conservation  of  energy  leads  to 
E(P  ,v  ,  X=l)+P  v  +1/2(D-U  )2  =  E  (P  ,l/x,A=0)+P  v  +1/2D2   (5) 

in  which  X  =  the  extent  of  chemical  reaction  (X  =0  for  no 
reaction  and  1  for  complete  reaction) .  Eliminating  u  and 
D  from  Eq  (5) ,  suing  Eqs  (1)  and  (2)  leads  to  the  Hugoniot 
curve : 
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B(P  ,v  ,X-1)  -   (Px,vx,X=0)  -  1/2  <Py+Px>  (Vx-v  )=0     (6) 

A  EH  =  1/2  (Py+Px) (vx-vy)  (7) 

AE    is  the  internal  energy  of  the  detonation  products  at  the 
very  high  pressure  and  temperature  of  the  detonation  state 
less  the  internal  energy  of  the  unreacted  explosive  at  ambient 
temperature  and  pressure.   It  may  be  written  as  the  sum  of  two 
terms : 

AEH  =  AET  +  AEp  (8) 

where   AE_   is  the  internal  energy  change  for  the  conversion 
of  reactants  to  products  isothermally ,  and   AE    is  that 
required  to  warm  and  compress  the  products  from  ambient   p 
and   T   to  the  detonation  values.   AE_   is  often  symbolized 
as  -Q  in  explosive  literature  (Q  being  the  heat  evolved  in  an 
explosion  when  corrected  to  isothermal  conditions) .   If  the 
products  were  ideal  gases,  then   AE   would  be  given   /c  dT, 
and  this  approximation  is  often  made,  though  the  products  at 
pressures  of  hundreds  of  bars,  are  far  from  being  ideal  gases. 

2.   The  Hugoniot  curve  is  that  of  a  hyperbola  in  the  p-v 
plane.   For  the  non-reacting  case  (X=0)  the  hyperbola  passes 
through  the  initial  state,  p  ,v  ;  but  for  the  case  we  are 
considering  here  (X=l)  it  does  not.   The  conservation  laws 
require  that  the  final  state,  however,  lies  on  the  Hugoniot 
curve  as  well  as  on  the  Rayleigh  line.   We  consider  below 
several  cases: 
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For  a  too  small  D(=d2)  the  Rayleigh  line  does  not 
intersect  the  Hugoniot  and  thus  no  solution  exists.   For 


Figure  2.8.   The  Hugoniot  Curve. 

large  D  (=D, ) ,  there  are  two  intersections,  called  strong 

(labelled  S)  and  weak  (W) .   For  one  particular  value  of 

D  (=D  .)  the  Rayleigh  line  is  tangent  to  the  Hugoniot. 

It  can  be  shown  that  the  relative  particle  velocity, 

D-u  ,  is  subsonic  at  S  (D=u   <  a  ) ,  supersonic  at   W  and 
P  P     Y "    * 

sonic  at  the  point  of  tangency,  called  the  Chapman- Jouguet 

point:   (D-u  )    =  a   (sonic  velocity  in  the  shocked  medium) . 

P  c]    y  * 

There  exist  a  number  of  computer  codes  based  on  such 
non-ideal  gas  treatments  (RUBY  and  TIGER  are  two  of  these)  and 
moderate  success  has  been  achieved  in  computing  detonation 
properties  from  these  codes.   The  equations  are: 
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Py   =   kpx2    <J>  (9) 

D      =    AcJ)*5    (1    +    BP    )  (10) 

<t>     =  nmV  (11) 


with   k  =  0.762,  A  =  22.3,  B  =  0.0013  (SI  units).   Q,  the 
heat  of  explosion  in  J/kg;  N,  the  number  of  moles  of  gas  per 
kg;  and  M,  the  average  molar  mass  of  the  gas  (kg/mole  are  all 
found  from  an  assumed  stoichiometry  of  decomposition,  combined 
with  thermochemical  principles.   Although  N,  M,  and  Q  are  all 
strongly  affected  by  the  assumption  made  as  to  stoichiometry, 
it  turns  out  that  the  important  parameter,  <J>,  is  rather 
insensitive  to  the  assumption. 
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III.   ANALYSIS  OF  DYNAMIC  RESPONSE  OF  PLATE  BY  SHOCK  WAVES 

A.   THEORETICAL  BACKGROUND 

1.   Nonlinear  Static  Analysis 

The  fundamental  method  used  to  obtain  a  static  solu- 
tion is  to  minimize  the  error  vector  {6},  given  by  the 
following  equation:   [Ref.  7] 

{6}  -  {P}+  {Q}-{F}  (12) 

where   {6}   is  the  error  vector  of  unbalanced  forces  acting 
at  all  grid  points,  {P}  is  the  vector  of  applied  external 
loads,  which  may  change  with  displacements,  {Q}  is  the  unknown 
vector  of  constraint  forces  due  to  single  and  multipoint  con- 
straints, and   {F}  is  the  vector  of  grid  point  forces  generated 
by  element  motion  and  stress.   The  terms  are  functions  of 
displacement,  temperature,  and  stress  history. 

It  must  be  noted  that  degrees  of  freedom  not  involved 
in  constraints  produce  null  terms  in   {Q}  and  that  dependent 
constraint  points  produce  no  errors.   Therefore,  when  Eq  (12) 
is  reduced  to  the  solution  coordinates  the  constraint  forces, 
{Q},  disappear.   Terms  in  the  vector   {F}  are  dependent  on  the 
deformations  of  the  finite  elements.   In  linear  static  analysis 
the  vector   {F}  becomes 

{Flinear}  =  ^^^  "  <PT(T)>  (13) 
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where  [1<P]    =  The  linear  stiffness  matrix 

T 
{P  }  =  The  "thermal  load"  vector,  i.e.,  a 

vector  of  grid  point  which  produces  the 

same  element  strains  as  the  temperature 

distribution, 

{u}   =  The  vector  of  grid  point  displacements. 

To  obtain  a  nonlinear  solution,  a  "tangent"  matrix 

[K ]  is  used  in  which  the  terms  are  derivatives  in  the  form 

6F. 
Kij  =  [*uHu  =  ur  ("> 

The  matrix   [K ]   contains  both  geometric  and  material  non- 
linear effects.   The  approximation  occurs  because  the  forces 
are  highly  variable  with  displacements  and  nonsymmetric  terms 
are  approximated.   If   {u  }   is  a  known  displacement  which 

produces  an  error   {5(u  )},  then  Newton's  method  may  be  used 

i+1 
to  predict  a  new  vector   {u    }   with  a  smaller  error.   For 

a  displacement   {u}   near  a  known  solution,   {u  },  the  non- 
linear force  is  approximately 

(F(u)}  *  (F(u1)}  +  [K]{u-u1}  (15) 

The  substitution  of  Eqn  (15)  into  Eqn  (12)  with   {6}  =  0 
yields  the  result 

[KHu-u1}  *  {P}  +  {Q}  -  (F(u1)}  (16) 

which  can  be  written  in  terms  of  the  error  {5}  to  provide 

the  format  of  the  Newton-Raphson  iteration  method: 

[K]{ui+1-ui>  =  (6(ui)}  (17) 
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When   [K]   may  be  inverted,  Eqn  (17)  may  be  used  to  solve 
for   {u    }  =  { (u) } ,  a  new  estimate  of  the  solution. 
Equation  (17)  provides  an  incremental  form  for  solution 
iteration.   Another  form  of  the  iteration  equation  is 
obtained  from  Eqns  (12)  and  (17) : 

[K]{ui+1-ui}  =  {pQ  +  pT}  +  {Q}  -  {f1}  -  [K]  {u1}    (18) 

where 

{f1}   =  {F(ui)}  -  [K]  {u1}  +  {pT}  -  Ap(ui)      (19) 

[K  ]  is  the  nonlinear  "reference"  matrix,  available  from  the 
elements  and   Ap(u1)   is  the  change  in  loads  due  to  grid 
point  motions.   The  advantage  to  this  form  is  that  the  second 
order  "corrective"  load  vector   {f  }  equals  zero  when  the 
forces   {F}   are  linear  and  when   [K]   contains  linear  elastic 
terms  (see  Eqn  13).   Because  the   {f}   forces  are  calculated 
at  the  element  level,  all  known  linear  elements  and  super- 
elements  may  be  bypassed  in  their  calculation.   This  is  very 
efficient  when  a  major  portion  of  the  structure  remains  linear 
However,  the  advantages  for  using  Eqn  (17)  are  that  the  load 
error  and  incremental  displacements   {u    -  u  }  are  directly 
available  for  testing  convergence.   Therefore,  in  MSC/NASTRAN, 
both  Eqns  (16)  and  (17)  are  used  internally. 

In  order  to  avoid  calculating  the  vector   {F(u  )}, 
Eqn  (19)  is  solved  for   {Ftu1)}   and  substituted  into  Eqn  (12) 
to  obtain 

{6i}  =  UCu1)}  =  {p  +  pT}  -   {f1}  -  [Kr]  {u1}       (20) 
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Subtracting  {6  }   for  two  successive  values  of   i   in 
Eqn  (20)  to  eliminate  the  constants  a  simple  equation  is  ob- 
tained for  the  load  error  in  the  form 

{61}  =  {5i"1}  -{f1  -  f1"1}  -  [k]  {u1  -  u1"1}    (21) 
This  form  provides  a  simple  iteration  procedure,  namely 

a)  For  the  first  iteration  after  a  change  in  the  load  vector, 
use  Eqn  (17)  and  (20) .   Note  that  the  initial  values  for 
starting  from  rest  are   {f}={u}=0:   calculate 

{5°}  =  {p}  +  {pT}   and  solve   [K]  {u1}  =   {5°}  (22)  (23) 

b)  then  compute    f    directly  from  the  element  and  load 
tables : 

{f1}   =  Flu1}  -  [Kr]  {u1}   +  {pT}  -  (Ap(u1)}    (24) 

c)  then  compute   {5  }   from  Eqn  (21) . 

{51}  =  6°  -  {f1}  +  {f°}  -   [Kr]  {u1  -  u°}       (25) 

where   {f  }  =  0  by  assumption. 

2 

d)  Use  Eqn  (17)  with  i=l  again  to  compute   {u  } 

[K]{u2  -  u1}  =  otu1}  (26) 

e)  Repeat  steps  b,  c  and  d  to  convergence,  incrementing  the 
superscripts  by  one  at  each  iteration  loop. 

2.   Nonlinear  Transient  Analysis  Using  MSC/NASTRAN 

SOL  99  which  is  a  rigid  format  of  MSC/NASTRAN  code 
for  nonlinear  transient  problem,  is  the  dynamic  comparison  to 
nonlinear  static  systems,  providing  for  combined  material  and 
geometric  nonlinear  analysis  in  the  time  domain.   SOL  66  which 
is  a  rigid  format  for  nonlinear  static  problem  and  SOL  99 
share  the  same  element  computer  code  and  solution  techniques 
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and  provide  similar  data  storage  and  restart  facilities. 
The  transient  solution  is  performed  in  a  stepwise  manner  with 
time  steps  replacing  load  steps  and  with  effects  of  mass, 
damping  and  unsymmetric  matrices  added  to  the  stiffness 
matrices . 

Because  the  calculations  are  identical,  many  element 
and  material  discussions  in  the  previous  sections  also  apply 
to  the  transient  analysis  system.   The  finite  elements  auto- 
matically provide  linear  inertial  forces  and  structural 
damping,  based  upon  the  initial  state,  as  well  as  nonlinear 
forces  due  to  geometry  and  material  effects. 

We  will  develop  the  basic  equations  for  the  transient 
solution.   The  methods  were  chosen  to  be  compatible  with 
existing  MSC/NASTRAN  static  nonlinear  solutions  as  well  as 
with  the  existing  linear  transient  methods.   Thus,  the 
Newmark  Beta  method  for  transient  integration  is  combined 
with  the  modifications  to  Newton's  method  for  nonlinear 
solutions.   The  additional  iteration  steps  provide  equilib- 
rium solutions  at  each  time  step,  thereby  guaranteeing 
stability  and  accuracy  for  arbitrary  time  step  size.   The 
basic  equations  are  developed  below,  followed  by  a  discus- 
sion of  the  stability  limits. 

a.   Basic  Equation 

Adding  dynamic  loads  to  the  basic  equation  given 
previously,  we  obtain  a  load  equilibrium  error  vector  {5  } 
at  time  step   n   by  equation: 
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{6}={P       -    Mu       -    Bu       -F 
IX  n  n  n  n 

where         ^pn^=    "Average"    load   over    the    time    period 

(t      ,     <  t       <  t_,) 

n-1  n  n+1 

•    •  • 

{U  }  {U  }=  Corresponding  acceleration  and  velocity  vectors, 
n    n  ■* 

{F  }=  "Average"  elasto-plastic  element  total  force 
n 

vector . 

The  above  equation  is  solved  at  the  "reduced" 

(u.)  displacement  vector  size.   The  approximation  errors  due 

to  dynamic  reduction  methods  are  not  included  in  the  error 

vector   {6  }.   Applying  Newmark ■ s  averaging  method  over  a 

finite  time  period,   t   ,  <  t  <  t  ,,   the  "static"  forces  are 
*  '        n=l        n+1 

{F  }  =   {$F(u  .-)  +  (l-2S)F(u  )  +  SF(u   ,)}         (28) 
n  n+l  n        n-l 

where   $  =  Newmark  Beta  operator  and   F(u  )   is  the  nonlinear 

force  due  to  a  generalized  displacement  vector  {U  }.   An 

identical  definition  occurs  for   {p  }   from  the  applied  loads 

^n  cc 

at  each  time  step. 

From  central  finite  differences,  the  acceleration 
and  velocity  vectors  are: 

{V  -  ^2  {iun+l  "  AV  -  ^2  {V  (29) 

(in>  "  2St  {Aun+l  +  iun}  (30) 

where    {Aun+1   =  {un+1  -  %} 

for  small  angles  and   At  =  t   ,  -  t    is  the  "time  step  size". 
3  n+1    n 

For  large  angle  changes  we  require  that   Cu} ,  {u},   and  {Au} 
are  vectors  in  the  global  coordinate  direction.   The  vector  {u} 

35 


contains  gimbal  angles  requiring  a  tranformation   R(u  ) 

such  that 

{u.}  =  {u  }   +  {R  }  Uu   .  }  (31) 

n+1       n       n      n+l 

the  load  error  in  terms  of  the  current  vectors  estimates: 

(V-<V  -  1-71 M  +  2itB)  {dn}  -1st-  biu%> 

At 
"   {Bfn+l  +  d-2S)Fn  +  trn_j)  (32) 


and 


{Au  .,  }  =  {Au   +  d  }  (33) 

n+l        n    n 


At  any  time  step  the  vectors  {P  },  {AU  }.  {F  }  ,  and  {F   ,  } 
1  r  n      n     n   '        n-1 

are  known.   In  general  the  vectors   {d  },  {u  .,}  and   f  ,,  = 

n     n+l        n+l 

F(u   ,)  must  be  found  either  by  approximating  {F   ,}   or  by 

an  iterative  search. 

Note  that  for  a  linear  solution  {F  }  equals 

n    ^ 

[K] (u  },  where  [K]  is  the  reduced  stiffness  matrix.   In  this 
n 

case  Eqns  (32)  and  (33)  may  be  used  to  solve  for  {u  ,,} 

*  n+l 

directly  by  setting  i&n}    =    0.   With   3=1/3  this  will  produce 
the  linear  SOL  27  MSC/NASTRAN  results. 
Nonlinear  Iterations 

For  a  nonlinear  solution,  Equns  (32)  and  (33)  may  be 
solved  with  Newton's  method  (or  a  modified  version) .   Using 
{d  }  as  the  primary  solution  variable  provides  the  following 
iteration  algorithm. 
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First  assume  a  linear  approximation 


t«aifi) 


{6n> 
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n 


3d 


n 


id 


i  +  1 


n 


O 


=  0 


(34) 


L-t 


{A1}  = 
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n 


3d 


n 


(35) 


Then,  the  estimated  displacement  change  is 


,i+l 


{d1+1  -  d1}  =  {A, }-i{5  (d1) } 
n      n      1     n  n 


(36) 


The  algorithm  is  identical  to  the  static  case  except  that 
instead  of  a  tangent  stiffness  matrix,  the  "left  hand  side" 
matrix,  obtained  from  Eqns   (32)  and  (34)  is 


{A,  }=  [-A*  M  +  to 


At 


2At 


(37) 


3F 


where   [K]  =  -*—  is  the  current  tangent  stiffness  matrix. 
At  a  new  time  step  the  load  iterations  may  be  started  by 

o 

assuming  that   {AU  . ,  }  =   {AU  }   and  therefore: 

n+1         n 


{6n>  ={0} 


(38) 


However,  to  be  consistent  with  displacements,  we  must  extra- 
polate the  nonlinear  forces: 


frtLlJ      =      {Fn     +      <Fn     -     F.      1>   }      =      ^2Fn     ~     Fn      1  }       (39) 

n+1      n     n    n-1        n    n-l 


Substituting  Eqns  (38)  and  (39)  into  Eqn  (32)  for  i=0  results 
in  the  first  estimate: 


^0=  (*J     +   [*,]   ^AUn>  "^Fn} 

n     n     L  2      n      n 


(40) 
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where 


[A2]  =   [-  jt  B]  (41) 


Here   A   is  a  precalculated  matrix.   {AU  }  and   {F  }   are 

n         n 

previous  results.   For  subsequent  iterations,   (i>0)  Eqn  (32) 
becomes : 

"n}  =  {0  -   1*3]^}  -  B(Fi+1-  F°+1}     (42) 

where  ,        , 

^  •    t^2  M  +  lit  BJ  (43) 

At  each  time  step  the  code  will  iterate  on  Eqns 

(36)  and  (42)  until    {6  }   passes  the  convergence  tests  or 

n   c 

the  number  of  passes  reaches  an  iteration  limit.   With  a 
single  step,  i=l,  calculating  only   6.,   the  results  will  be 
identical  to  the  existing  NOLIN  results  in  MSC/NASTRAN.   For 
faster  convergence  the  iterations  may  continue,  the  matrices 
may  be  updated,  and/or  the  time  step  size  may  be  reduced. 
3 .   Factors  Determining  Response  to  Transient  Loading 
The  response  of  a  metal  to  dynamic  loading  depends 
upon  three  elements  which  are  characteristic  of  the  material: 
the  dynamic  constitutive  relation,  the  pressure-temperature 
phase  diagram,  and  the  dynamic  fracture  criteria.   In  the 
absence  of  phase  changes  or  fracture,  the  material  response 
is  governed  completely  by  the  constitutive  relation,  or 
dynamic  equation  of  state.   The  constitutive  relation  is 
characterized  as  "elastic",  "viscoelastic" ,  "visco-elastic- 
plastic" ,  etc.   It  depends  on  macroscopic  parameters  such  as 
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Young's  modulus  and  the  viscosity  coefficient  which,  on  a 
microscopic  scale,  arise  from  atomic  interactions  and  dislo- 
cation processes.   The  phase  diagram  determines  whether  a 
given  loading  history  will  bring  the  metal  into  a  pressure- 
temperature  regime  where  phase  transitions  are  possible.   If 
this  occurs  during  shock  loading,  multiple  shock  waves  may 
result.   Microcracks  and/or  voids  begin  to  appear  in  a  metal 
under  tension  whenever  the  criteria  for  nucleation  and  growth 
of  damage  are  met.   These  criteria  are  functions  of  material 
parameters  and  the  stress  history. 

In  this  chapter  the  dynamic  constitutive  relations 
are  discussed.   It  is  the  most  important  element  to  determine 
the  response  of  a  metal  to  dynamic  loadings. 

In  macroscopic  or  rheological  language,  iron  can  be 
described  as  a  viscoelastic-plastic  material,  i.e.,  it 
responds  elastically  until  yielding  (with  upper  and  lower 
yield  points) ,  and  is  strongly  rate-dependent. 

The  constitutive  relation,  or  dynamic  equation  of 
state,  thus  gives  the  stress  as  a  function  of  strain  and 
time,  or,  alternatively,  can  be  written  as  a  differential 
equation  in  o,z,o,£.      and  an  important  goal  of  present  day 
research  is  to  develop  the  capability  to  derive  this  macro- 
scopic constitutive  relation  by  averaging  the  effects  of 
microscopic  processes. 

This  is  analogous  to  the  way  in  which  the  thermo- 
dynamic equation  of  state  of  a  gas  is  derived  by  averaging 
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the  effects  of  the  molecular  collisions.   Once  developed, 
this  capability  would  be  an  extremely  powerful  tool,  since 
the  response  of  a  metal  to  dynamic  loads  could  be  predicted 
from  known  microscopic  material  properties.   Although  this 
goal  has  not  been  attained  as  yet,  there  has  been  significant 
progress  during  the  last  few  years  in  this  direction  for  the 
case  of  crystalline  materials,  in  which  dislocation  processes 

play  the  role  of  the  "molecular  collisions".   In  iron,  the 

8    2 
density  of  mobile  dislocations  is  of  the  order  of  10  /cm  , 

so  the  number  of  dislocations  present  is  certainly  large 

enough  to  make  a  statistical  approach  valid. 

The  basic  approach  is  based  on  the  relation  from 

dislocation  theory:  [Ref.  3] 

YP  =  N   bv  (44) 

m 

where   y   is  the  plastic  shear  strain  rate 

(y-  •  =  t  (^ +  ^"^  )  /   where   £.   is  the  displacement  in 

11     2.         dX  .     dX .  1 

the  x.  direction) ,   N    is  the  mobile  disclocation  density, 
b  is  the  Burgers  vector,  and   v   is  the  average  disclocation 
velocity.   This  equation  provides  the  link  between  microscopic 
dislocation  motion  and  macroscopic  plastic  straint  rate,  and 

is  based  on  the  fact  that  it  has  been  found  possible  for  many 

p 

crystalline  solids  to  express   N   and   v   as  functions  of  y 

■*  m 

and  applied  macrosopic  shear  stress,   t. 
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The  constitutive  relation  is  obtained  by  combining 
Eqn  (44)  with  the  incremental  Hooke ' s  law 

da.  .  =  2ydee.  +  AdE®   6.  .  (45) 

under  the  assumption  that  the  total  strain  is  the  sum  of 
elastic  and  plastic  components: 

ij    £ij  +   £ij  (46) 

where  A   and    u   are  the  Lame  constants.   5. .   is  the 
Kronecker  delta,  and  a  repeated  index  indicates  summation. 
The  plastic  strain  is  assumed  to  cause  no  volume  change,  i.e., 

£KK  "  °  (47» 

The  combination  of  Eqns  (46) -(48)  yields,  after  taking  the 

time  derivative: 

a.  .  -2  we.  .  -  Aevv5.  .  -  2  ye?.  (43) 

ij  ID      KK  i]        13 

In  order  to  combine  Eqn  (48)  with  Eqn  (44) ,  we  must 
first  specify  the  loading  geometry.   In  plate-slap  tests 
under  conditions  of  uniaxial  strain,  all  total  strain  com- 
ponents except  e    in  the  direction  of  wave  propagation  are 
zero,  and   £-— -  e  .   Then  equation  (48)  becomes 

a   -  (A+2u) e   =  -2  y  eP  (49) 

XX  X 

Also,  in  this  case,  the  plane  of  maximum  shear  stress  in 
isotropic  media  is  inclined  at  45°  to  the  direction  of  shock 
propagation,  and 
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The  dynamic  constitutive  relation  for  this  loading  geometry 
is  thus  obtained  by  combining  Eqn  (44) ,     (49) ,  and  (50) : 

or   -    (A    +    2M)£V    =    -   4  uNbv  (51) 

Vv  J\  O  ill 

If,  A,   u,   and  the  right-hand  side  of  Eqn  (51)  are  known 
as  functions  of  macroscopic  stress  and  strain,  then  the 
equation  can,  in  principle,  be  solved  for  any  problem 
involving  uniaxial  strain. 

4 .   Dynamic  Plastic  Deformation  of  Plates 

We  consider  cases  of  the  impulsive  loading  of  thin 
metal  plates,  a  simple  situation  in  which  an  exponential 
shock  wave  impinges  on  a  clamped  square  plate.   In  this  case, 
the  plastic  deformation  is  supposed  to  be  brought  about  by 
bending  only  and  combined  bending  and  tension  is  not  analy- 
tically considered.   The  pressure  in  the  shock  wave  may  be 
assumed  to  be  given  by  [Ref.  4], 

P  =  Pm  exp(-t/x) 
where   p   is  the  peak  pressure  at  a  point  at  the  head  of  the 

wave,  t  denotes  time  and   x   is  a  time  constant  the  time 

taken  for  the  pressure  to  fall  to   p  ,e ,    which  for  high 
explosives  is  ~30   usees.   We  obtain  simply,  some  of  the 
results  of  Cox  and  Morland  [Ref.  3  ]f  who  theoretically 
investigated  the  dynamic  deformation  of  plates,  side  length 
2a,  position-fixed  at  their  periphery,  when  subject  to  a 
uniform  constant  transverse  pressure   p   for  time   x.   Their 
detailed  examination  showed  that  two  categories  of  load  may 
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be  distinguished;  medium  load  when   p  <   p  <  2p    and  high 
load  when   2p  <   p,  where   p    denotes  the  pressure  required 
just  to  enforce  plastic  collapse  of  the  plate,  statically 
For   p   we  have,  see  Fig.  3-1. 


4a/2~  •  W/2  .  M  =  i  •  4a2  -aw  •  P 

P    3  s 


(52) 


i.e. 


P   =  6Mn/a' 
s      P' 


For  simplicity  we  deal  only  with  their  first  case, 
and  for  the  second,  more  complicated  one,  the  reader  is 
referred  to  the  original  source.   For  this  medium  range  of 
loading  two  sub-divisions  are  required:   (1)  when   0<t<i 
and  (2)  when  t>x.    Use  is  made  of  the  results  obtained 
above  for  static  situations  using  hinge  lines,  and  we  take 
over  whatever  results  may  be  useful.   As  previously,  the 
analysis  pays  no  special  attention  to  material  rate  effects 
Further,  it  is  implicit  in  this  approach  that  the  static 
and  dynamic  modes  of  deformation  are  identical. 


NET  ON  THE 
FRAME 


a 
< 


<d 


TIME 


TIME 


Figure  3.1.   Dynamic  Loading  of  Square  Plate 
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a.   Medium  Load:   0  <  t  <  T 

Assume  the  mode  of  deformation  of  the  plate  is 
that  of  Fig.  3.2  and  that  the  plate  diagonals  become  hinge 
lines.   Since   p   is  constant  and  greater  than   p  ,   and 
because  the  four  equal  parts  of  the  plate  act  as  rigid  bodies 
each  rotating  about  its  boundary  or  side,  they  must  undergo 
an  angular  acceleration,  a.   If  changes  in  plate  geometry 
throughout  the  deformation  are  neglected,  the  resistance  to 
deformation,  which  arises  from  the  yielding  taking  place  in 
the  hinge  lines,  will  be  constant;  this  resistance  is  related 
to   p  ,  the  pressure  to  cause  static  yield.   Thus  the  excess 
load  is  constant  and  hence   a   must  be  constant.   The  work 
done  by  the  excess  load  will  manifest  itself  as  kinetic 

energy  acquired  by  the  four  rigied  regions  of  the  plate,  i.e. 

1    ■> 
4  *  y  Ioj  ,   where   I  is  the  moment  of  inertia  of  a  triangle 

about  a  side  of  the  plate  and   co   is  its  current  angular 

velocity.   Because   a   is  constant,  at  the  end  of  time  the 

plate  centre  will  have  descended  through  a  distance   vt/2, 

where   v   is  its  current  speed.   Thus,  an  equation  involving 

the  work  available  for  giving  rise  to  rotational  kinetic 

energy  of  the  plate  triangles  is, 

2    vt/2        1 
(p  -  p  )  4az  * =4  •  i  Ioj2  (53) 

S  3  Z 

4 
Hence  substituting  for   p    from  (52) ,  noting  that  I  =  ma  /6 

3 


and  simplifying, 
I 
2  m 
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6M                      (p-p   t 
(p  -  — ^— )  t  =    /2  or  v  =  2  - (54) 

a 


or  directly  by  momentum  consideration,  thus, 

4 
s  6     a 


,  ,     ,        a •  2a,  a    ma 

((p-p  )  •  t  •  -* )T  =  — =—  •  f  (55) 


The  deflection   w,   at  the  end  of  time   x   is, 

wx  =  VT/2  =  (p  -  ps)i2/m  (56) 

(2)  Medium  Load:   t  >  t 

After  the  removal  of  pressure   p   at   t  =  x,  further  deflec- 
tion,  w_ ,   is  added  to   w,,   because  the  rotational  kinetic 
energy  of  the  four  plate  segments  must  be  dissipated  in  doing 
plastic  work  in  the  hinges.   Let  the  segments  rotate  through 
further  angle  $      before  coming  to  rest,  then, 

1     v2  w?/ 

4  •  (^  I  * — T)  =  4  a/2  •  M   (<J>/2)  =  8aM   ■  ^/ a 

a2  p  p 

Hence 

2Tr2        2    (p-p  )  2t2 
ma  v       mv   _   ^  rs  ,_... 

W2    2  4  M   "   4p   "     m  p  [b  '' 

p     rs        ^s 

using  Eqn  (55).   Thus  the  total  central  deflection  is, 

2 

w  =  w1  +  w2  =  p  (p  -  pg) x  /mps  (58) 

If  the  total  response  time  of  the  plate  is   t  ,   note  that 
x.p   =  px.   In  the  case  of  thin  plates  clamped  at  their  outer 
periphery  it  is  clear  that  after  a  significant  amount  of 
deflection  there  must  be  stretching.   An  analysis  of  the 
dishing  of  peripherally  clamped  metal  plates  when  given  a 
transverse  initial  speed  of  a  few  hundred  feet  per  second 
has  been  given  by  Hudson  [Ref.   3].   Descriptively,  Hudson's 
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model  is  easily  understood  by  reference  to  Fig.  3-2.   Central 
portions  of  the  plate  maintain  their  given  transverse  speed 
until  overtaken  by  a  bending  wave  propagated  radially  inwards 
from  the  clamping  ring.   Fig.  3-2 (a)  shows  the  imagined  plate 
configuration  at  successive  intervals  during  the  deformation 
process;  when  the  plate  deflection  is  significant,  the  process 
is  predominantly  one  of  stretching  and  for  many  practical 
purposes  the  technological  approach  is  adequate. 

Figure  3-2 (a)  shows  the  terminal  state  taken  up 
by  a  thin  rectangular  steel  plate  clamped  around  its  periphery 
after  being  subjected  to  the  uniform  but  intense  impulsive 
loading,  the  maximum  speed  of  the  plate  is  a  few  hundred 
feet  per  second,  but  contracting  in  area  with  time;  the 
dynamic  situation  during  deformation,  with  plastic  hinges 
moving  into  the  plate  from  the  clamped  edges  and  intersecting 
along  lines  equally  inclined  to  the  sides  from  the  corners, 
is  indicated  diagrammatically  in  Fig.  3-2 (b). 

Now  the  yield  criteria  of  metals  should  be  con- 
sidered.  Two  approaches  to  the  interpretation  of  a  yield 
criterion  are  possible;  one  is  Tresca's  criterion,  the  other 
is  Von  Wises  criterion.   But  for  most  metals  it  is  generally 
found  that  while  the  experimental  points  fall  between  the  two 
ellipses,  they  incline  towards  Von  Mises '  ellipse.   Experiments 
were  earlier  performed  by  Lode  (1925)  on  tubes  subject  to 
tension  and  internal  pressure.   Bending  and  torsion  was  used 
by  Siebel  in  1953.   All  results  confirm  the  general  opinion 
that  the  best  simple  yield  criterion  for  metals  is  that  of 

Von  Mises  [Ref.  51.  .  _ 
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Figure  3.2.   Dynamic  Deformation  of  Fixed  Square  Plate 
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B.   PROBLEM  DESCRIPTION  AND  NASTRAN  DECK  SETUP 
1 .   Problem  Description 

The  model  is  shown  in  Fig.  3.3.   It  is  a  20  in.  by 
20  in.  steel  square  plate  model  which  is  divided  to  25  QUAD4 
(NASTRAN  card)  elements  and  its  boundary  is  clamped.   Its 
thickness  is  0.8  inches  and  the  engineering  data  as  follows. 

E  =  3.0E+7 

v  =  0.3 

L  =  20.0  in. 

p  =  7.33E-4  lb/in3 

t  =  0.8  in. 
where   E   is  the  elastic  modulus,   v   is  the  Poisson's  ratio, 
L  is  the  length  of  the  plate,   p   is  the  mass  density,   and 
t  is  a  thickness  of  the  plate. 

Next  the  shock  wave  applied  impinges  on  the  steel 
plate  in  the  -Z  direction  at  time  0.0  seconds  (Fig.  3.4)  and 
its  peak  overpressure  is  3.0E+4  psi  and  the  shock  wave  dura- 
tion time  (t,)  is  5.0E-3  sees.   To  find  the  pressure,  the 
following  formula  is  used  [Ref.  6] 

P  =  Pm^-t/V  exp(-at/td) 

where   p    is  peak  reflected  overpressure  at  time  0.0  sees 
and   a   is  wave  form  parameter.   But  the  shock  wave  pressure 
cannot  be  applied  as  a  load  to  the  plate  since  the  front  face 
shock  pressure  of  the  plate  is  relieved  by  the  rear  face  shock 
wave  which  is  positive  Z  direction.   The  following  formula  is 
used  [Ref.  4  ]  . 
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p  =  p   exp(-t/x) 

r     m 

The  load  data  as  follows. 

peak  lateral  loads  (-Z  direction) 
p  =  3  .E+4  psi  on  the  Element  (213) 

p  =  2.2E+4  psi  on  the  Element  (207,  208,  209,  212, 
214,  217,  218,  219) 
Von  Mises  criterion  would  be  applied  for  this  model.   The 
time  step  (t)  to  be  applied  is  1.0E-6  sees  from  0.0  sec  to 
7.0E-4  sees. 
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Figure  3.3.   Modeled  Steel  Square  Plate. 
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Figure  3.4   Shock  Wave  Applied. 
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2.   NASTRAN  Input  Data  Deck 

Reference  is  made  to  the  echo  of  the  input  data  deck 
which  appears  in  the  output, 
a.   Sol  66  [Ref.  6] 

1  NASTRAN  PREF0PT=2 

2  ID  NONLINEAR,  STATIC 

3  SOL  66 

4  TIME  3  0 

5  CEND 

6  TITLE-NONLINEAR  DEFORMATION  WITH  STATIC  LOADS 

7  SUBTITLE=TRANSIENT  RESPONSE  BY  SHOCK  WAVE  LOAD. 

8  ECHO=UNSORT 

10  DISP=ALL 

11  OLOAD=ALL 

12  SPCF=ALL 

13  GPFOR=ALL 

14  ELST=ALL 

15  SEALL=ALL 

16  SPC=200 

17  SUBCASE  1 

18  LOAD=10 

19  NLPARM=10 

20  BEGIN  BULK 

21  $ 

22  $  LATERALLY  LOADED  FLAT  PLATE 

23  $ 
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24  MESHOPT, , , , ,YES 

25  $ 

26  $  USE  MSGMESH  TO  GENERATE  5*5  ARRAY  OF  QUAD4  ELEMENTS 

27  $ 

28  EGRID,1, ,0. ,0. 

29  EGRID, 2, ,20. ,0. 

30  EGRID, 3, ,20. , 20. 

31  EGRID, 4, ,0. ,20. 

32  GRIDG, 2, , ,5,-1,-2,-3, , +FLD2 

33  +FLD2,5,-4 

34  CGEN,QUAD4,201,200,2 

35  $  REFER  TO  PAGE  2 . 18-7, APPLICATION  MANUAL. 

36  PSHELL,200,65,0.8,  65, ,65 

37  MATl,65,3.0E+7, , . 3,7.33E-4 

38  MATS1, 65, 40, PLASTIC, ,1,1,2. 6E+4 

39  TABLES  1 , 40 , 0 ,,,,,,, +TAB1 

40  +TAB1 , 0.0,0.0,0.001,3. 0E+4 ,0.025,4. -E+4 ,0.05,5. 8E+4 , +TAB2 

41  +TAB2 ,0.10,6. 7E+4 ,0.125,7. OE+4 ,0.2,7. 5E+4 , ENDT 

42  PARAM,COUPMASS,l 

43  PARAM,K6ROT,1.0 

44  PARAM,LGDISP,1 

45  PLOAD2, 10, -14. 7, 2 0  8, 2 12, 2 13, 214, 218 

46  NLPARM,10,1, ,AUTO, ,20,W,YES 

47  $NLPARM,20,3, ,AUTO, ,20,W,YES 

48  SPCG,200,2,123456,A,B 

49  SPCG, 200, 2, 123456, B,C 
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50  SPCG,200,2, 123456,0,0 

51  SPCG, 200,2, 123456, A, D 

52  ENDDATA 

1)  Executive  control  card 

card  1  -  The  NASTRAN  card  contains  the  keyword   PREFOPT 
=  2  which  requests  that  MSGMESH  be  used  to 
generate  input  data. 

card  3  -  Request  nonlinear  static  analysis,  Rigid  format 
66. 

2)  Case  control  card 

card  15  -  Requests  superelement  generation  and  assembly 

*  should  be  set  as  'all' 

card  18  -  Requests  static  load  combination. 

*  SID  should  be  matched  to  the  PL0AD2  card  in 
the  bulk  data  deck. 

card  19  -  Requests  nonlinear  parameter  selection. 

*  SID  should  be  matched  to  the  NLPARM  card 
in  the  bulk  data  deck. 

*  This  card  may  appear  above  or  within  a 
subcase. 

*  LOAD  card  is  only  applicable  in  statics, 
buckling,  heat  transfer  problems. 

3)  Bulk  data  cards 

card  24  -  Requests  out  options  control 

*  An  entry  of  YES  causes  the  printing  of  the 

*  input  card  echo ' . 
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card  28-31  -  Defines  vertices  of  quadrilateral  field 

using  MSGMESH  EGRID  command, 
card  33  -   Defines  a  field  of  grid  points  and  generate 
grid  cards. 

*  The  grid  number  generated  by  MSGMESH  starts 
from  the  FID  of  GRIDG  (if  FID  is  2,  the  first 
grid  point  ID  is  20000) 

*  5  means  5  elements  in  each  direction, 
card  34  -  Generate  element  connection. 

*  QUAD4  is  the  element  type  which  should  be 
generated. 

*  Element  ID  number  starts  from  201. 

*  Fifth  field  200  is  the  PID  of  pshell  card, 
card  38  -  Specifies  table  references  and  material 

properties  which  are  stress-dependent  for  use 
in  material  nonlinearity  applications 

*  Taken  Von  Mises  yield  function,  isotropic 
hardening  rule. 

card  39  -  Represent  stress-strain  data. 

card  43  -  The  value  1.0  suppresses  singularities 

*  It  is  intended  primarily  for  geometric 
nonlinear  analysis. 

card  44  -  The  value  1  allows  the  nonlinear  calculations, 
card  45  -  Request  the  initial  condition  on  the  elements 

in  the  negative  Z  direction, 
card  46  -  Controls  the  nonlinear  analysis  iteration. 
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*  To  reduce  the  output,  use  sixth  field. 

*  For  rapid  convergence,  increase  the  Error 
tolerance  and  maximum  iteration  and  the  limit 
on  diverging  iterations. 

card  4  8-51  -  Generate  spcl  cards. 

This  card  is  only  used  to  MSGMESH  generation, 
b.   Sol  99 

1  NASTRAN  PREF0PT=2 

2  ID  NONLINEAR, TRANSIENT 

3  SOL  99 

4  TIME  59 

5  CEND 

6  TITLE=NONLINEAR  DEFORMATION  WITH  DYNAMIC  LOADS 

7  SUBTITLE=TRANSIENT  RESPONSE  BY  SHOCK  WAVE  LOAD. 

8  ECHO=SORT 

9  SET  1=20202,20203,20204,20304,20404 

10  SET  2=201,208,212,213,214,219,225 

11  SET  3=20000,20003,20005,20200,20400 

12  SDISP=1 

13  OLOAD=l 

14  SPCF=3 

15  ELST=2 

16  SEALL=ALL 

17  SPC=200 

18  LOADSET=77 

19  SUBCASE  1 


55 


20  DLOAD=3  7 

21  TSTEPNL=45 

22  SUBCASE    2 
2  3  DLOAD=3  7 

24  TSTEPNL=46 

25  SUBCASE  3 
2  6  DLOAD=3  7 
27  TSTEPNL=47 

2  8  OUTPUT (XYPLOT) 

29  PLOTTER  NAST 

30  CSCALE  1.3 

31  XTITLE=TIME  IN  SECS 

32  XGRID  LINES=YES 

33  YGRID  LINES=YES 

34  YTITLE=DISPLACEiMENT  GRID  20202 

35  XYPLOT  SDISP  RESP/20202 (T3) 

36  YTITLE=VELOCITY  GRID  20202 

37  XYPLOT  SVELO  RESP/20202 (T3) 

38  YTITLE=DISPLACEMENT  GRID  20101  20102  20202 

39  XYPLOT  SDISP  RESP/20101 (T3) , 20102 (T3) , 20202 (T3) 

40  BEGIN  BULK 

41  $ 

42  $  LATERALLY  LOADED  FLAT  PLATE 

43  $ 

44  MESHOPT, , , , ,YES 

45  $ 
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4  6  $  USE  MSGMESH  TO  GENERATE  5*5  ARRAY  OF  QUAD4  ELEMENTS 

47  $ 

48  EGRID, 1, ,0. ,0. 

49  EGRID, 2, ,20.  ,0. 

50  EGRID, 3, ,20. ,20 . 

51  EGRID, 4, ,0.  ,20. 

52  GRIDG, 2, , ,5,-1,-2,-3, , +FLD2 

53  +FLD2,5,-4 

54  CGEN,QUAD4,201,200,2 

55  $    REFER   TO    PAGE    2 . 18-7  APPLICATION    MANUAL 

56  PSHELL,200,65,0.8,65,  ,65 

5  7  MATl,65,3.0E+7, , .3,7.3  3E-4, , ,0.01 

58  MATS1, 65, 40, PLASTIC, ,1,1,2. 6E+4 

59  TABLES1,40,0, , , , , , ,+TABl 

60  +TAB1 ,0.0,0.0,0.001,3. OE+4 ,0.025,4. OE+4 ,0.05,5. 8E+4 , +TAB2 

61  +TAB2,0.10,6 . 7E  +  4 , 0 . 125 , 7 . OE+4 , 0 . 2 , 7 . 5E+4 , ENDT 

62  PARAM,COUPMASS,l 

6  3  $PARAM, AUTOS PC, YES 

64  PARAM, W4,1.616E+3 

65  PARAM,K6ROT,1.0 

66  PARAM,  LGDISP,1 
6  7  PARAM,  SLOOPID,l 
6  8  PARAM,LOOPID,l 

69  PARAM,  STIME, 0.000450 

70  PARAM,MAXLP,7 

7 1  DLOAD ,37,1.0,1.0,7,1.0,6,1.0,5 
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72  TL0AD1,7,61,  ,0,22 

73  LSEQ, 77,61,71 

74  PLOAD2,71,3.2E+4,213 

75  TL0AD1,6,62, ,0,22 

76  LSEQ, 77, 62, 81 

7  7  PLOAD,81,2.5E+4,20  7,20  8,209/212,214 

78  TLO AD1 , 5 , 6 3 ,  ,0,22 

79  LSEQ, 77, 63,91 

80  PLOAD2, 91, 2. 5E+4f217, 218,219 

81  TABLED1, 22, ,,,,,, , +NEXT1 

82  +NEXT1,-1. ,0.0,0.0,-0.0  010,0.0  001,-1.0,0.0  010,-0.56  0,+NEXT2 

83  +NEXT2, 0.004,-0. 001, ENDT 

84  TSTEPNL,4  5,15  0,1.0E-6,20,AUTO, , 20,W,+LLl 

85  +LLl,1.0E-7,1.0E-7, 1.0E-7,20,20,10 

86  TSTEPNL,46,200,1.0E-6,20,AUTO, ,20,W,+LL2 

8  7  +LL2,1.0E-5,1.0E-5,1.0E-5,20,20,10 

88  TSTEPNL,47/120,3.0E-6,20,AUTO/ ,20,W,+LL2 

89  +LL2,1.0E-4,1.0E-4,  1. OE-4 , 20 , 20 , 10 

90  SPCG, 200,2, 123456, A, B 

91  SPCG, 200, 2,123456, B,C 

92  SPCG, 200, 2, 123456, C,D 

93  SPCG, 200, 2, 123456, A, D 

94  END DATA 

1)  Executive  control  cards 

card  3  :   Requests  nonlinear  transient  analysis,  Rigid 
format  99. 
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card  25-27 
card  28-39 


2)  Case  control  cards 

card  18:   Selects  a  static  load  set  for  use  in  dynamics. 

Thus  this  may  be  referenced  by  dynamic  loads, 
card  20 :   This  card  selects  a  dynamic  load  for  nonlinear 

transient  problem, 
card  22-24:   Requests  restart  at  time  0.00015  sees. 

Requests  restart  at  time  0.00035  sees. 
Requests  the  xy  plotting  which  is 
corresponding  to  the  time  step. 

3)  Bulk  data  cards 

card  64:   This  gives  damping  effects  for  transient 
analysis . 

*  Sometimes,  more  rapid  convergence  effects 
can  be  obtained  by  this  damping  factor  with 
damping  coefficient  in  MAT1  card  (GE) . 

card  67:   This  card  identifies  the  initial  conditions 
from  a  previous  SOL  66  nonlinear  static 
solution. 

*  Refer  to  the  param  card,  SLOOPID. 

card  6  8,69:   Request  these  two  cards  when  we  taken  more 

than  one  subcase . 
card  70:   Requests  more  internal  calculation. 

*  We  can  reduce  iteration  steps  and  prevents 
the  overprinting  due  to  iteration  steps. 
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card  71:  Defines  a  dynamic  loading  condition  for 
transient  response  problems  as  a  linear 
combination  of  load  sets 

*  See  the  remarks  5,  7,  9  of  DLOAD  card  in 
bulk  data. 

card  72,  75,  78:   Defines  time  dependent  dynamic  load. 

*  Requests  a  static  load  card  set  (LESQ) . 
card  73,  76,  79:   Defines  a  sequence  of  static  load  sets 

This  load  sets  may  be  referenced  by  dynamic 
load  card  (TL0AD1) 

*  This  card  should  be  used  when  we  want  to 
apply  area  pressure  on  elements  than  a  point 
pressure  on  grid  point. 

card  74,  77,  80:   Defines  a  uniform  static  pressure  load 

on  the  plate  elements, 
card  81-83:   Defines  a  tabular  function  for  use  in 

generating  a  time-dependent  dynamic  loads. 

*  In  this  card  the  peak  pressure  (3.0E+4  psi) 
occurred  at  0.0001  sec. 

card  84,  86,  88:   This  card  control  the  nonlinear 
transient  analysis. 

*  This  card  is  expalined  in  the  sensitivity 
analysis. 

C.   SENSITIVITY  ANALYSIS 

1.   General 

In  transient  analysis,  two  types  of  instability  could 

occur.   The  first  is  the  familiar  nonlinear  load  iteration 

divergence  which  also  occurs  in  static  analysis.   The  second 
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is  the  divergence  which  grows  with  time  in  the  transient 
integration.   Both  instabilities  are  caused  by  uncorrected 
nonlinear  equilibrium  errors. 

A  convenient  method,  from  Von  Neumann,  for  analyzing 
the  stability  limits  is  to  assume  that  the  nonlinear  forces 
are  nearly  linear  and  the  error  vector  has  a  constant  conver- 
gence.  It  is  assumed  that  the  nonlinearity  has  a  first  order 
approximation:  [Ref.  7] 

(F(y) }  -  [Kr  +  AKNL]  {y}  (59) 

NL 
where   AK     is  the  difference  between  the  tangent  stiffness 

matrix  and  its  approximation   K  .   The  error  vectors,   {6}, 

are  assumed  to  grow  at  the  rate   A,  defined  as: 

{5i+1}=  A{6.}  (60) 

Note  that  if   |X|  >  1   the  system  will  be  defined  as  unstable. 

After  lengthy  calculations  we  may  summarize  the  various 
criteria,  assuming  that  the  matrices  are  reduced  to  equivalent 
scalar  modal  quantities,  the  criterian  for  stable  solutions 
are : 

1)  For  time  steps  with  converged  "static"  iterations: 

1/2  >  S  >  1/4 
This  is  the  same  criteria  as  linear  analysis. 

2)  For  "static  load"  iterations: 

3  A  KNL  <  8  Kr  +  m  B  +  -ij-  M 
~        ZAC     At 

This  states  that  the  mass  and  damping  add  to  the  effective 
linear  stiffness  [Kr]  and  improve  the  stability. 
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3)  For  time  step  integration  with  no  intermediate  static 
iterations,  as  in  the  standard  Newmark  Beta  method: 

NL  r     4 

A  K    <_    (4  0  -1)  Kr  +  — *-*■   M 

At 
This  restriction  is  more  severe  than  criteria  b  above, 
proving  that  the  internal  iterations  are  more  stable 
than  the  Newmark  integration. 

In  summary,  the  Von  Neumann  method  will  have  fewer 
divergence  problems  than  either  the  static  nonlinear  solu- 
tion or  the  single  step  transient  nonlinear  methods.   The 
better  stability  of  the  method  and  the  capability  to  use 
larger  time  steps  outweighs  the  cost  of  a  few  internal 
iterations  on  the  static  element  forces. 
2 .   Nonlinear  Transient  Controls 

The  input  data  required  for  SOL  99  is  a  combination 
of  transient  control  data,  similar  to  SOL  69  (Direct  Transient, 
superelements) ,  and  nonlinear  modeling  data  similar  to  SOL  66 
(Nonlinear  Statics) .   The  nonlinear  properties  are  defined  by 
nonlinear  material  data  (MATS1) ,  gap  elements  (GAP) ,  and 
large  displacement  effects  (parameter  LGDISP) .   The  transient 
effects  are  produced  by  loading  functions  (TLOADi,  DAREA,  etc.), 
damping  (parameters,  elements,  and  material  data),  and  element 
mass  properties. 

The  only  unique  data  required  for  SOL  99  is  supplied 
on  the  TSTEPNL  data  cards.   (Several  optional  parameters  are 
provided  for  tuning  the  method  and  for  use  in  restarts.)   The 
TSTEPNL  card  in  itself  is  an  amalgam  of  the  normal  transient 

TSTEP  card  and  the  nonlinear  control  card,  NLPARM. 
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a.   Case  Control 

The  significant  rules  for  nonlinear  transient  are: 

1)  Each  subcase  defines  a  time  interval  starting  from  the 
last  time  step  of  the  previous  subcase  which  is  subdivided 
into  small  time  steps.   The  output  TIME  printout  will  be 
labeled  by  the  net  time,  including  all  previous  subcases. 

2)  The  input  loading  functions  may  be  changed  for  each  sub- 
case or  continued  by  repeating  the  same  DLOAD  request. 
The  bulk  functions  are  evaluated  at  the  net  model  time, 
which  starts  at  the  last  time  step  of  the  previous  subcase. 

3)  Each  subcase  may  have  a  different  time  step  size,  time 
interval,  and  iteration  control  selected  by  the  TSTEPNL 
request.   The  case  control  requests  which  may  not  be 
changed  after  the  first  subcase  are:   SPC,  MPC  ,  DMIG,  and 
TF. 

4)  Output  requests  for  each  subcase  are  processed  indepen- 
dently and  the  output  for  each  case  are  generated  before 
the  next  subcase  is  processed.   Separate  XYPLOT  ouputs  will 
occur  for  each  subcase.   Available  outputs  are  DISPLACEMENT, 
VELOCITY,  ACCELERATION,  OLOAD ,  STRESS,  FORCE,  S DISPLACEMENT, 
SVELOCITY,  and,  SACCELERATION .   Note  that  NONLINEAR  and 
SPCFORCE  outputs  are  not  provided. 

5)  Initial  conditions  are  not  explicitly  provided  for  user 
input  in  SOL  99.   The  reason  is  that  the  nonlinear  element 
configuration  and  equilibrium  state  must  be  known  at  every 
time  step.   If  initial  conditions  were  given,  all  of  the 
nonlinear  element  forces  and  plastic  stresses,  as  well  as 
the  past  history,  must  be  given.   Therefore,  the  initial 
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conditions  become  part  of  the  analysis  and  are  generated 
by  initial  transient  subcase  or  by  restarts  from  a  previous 
static  solution. 
6)  All  normal  superelement  model  generation  options  and 

matrix  reduction  options  are  allowed  for  the  linear  portion 
of  the  structure.   General  Dynamic  Reduction,  Component 
Mode  Synthesis,  and  Guyan  Reduction  may  be  performed  for 
nonresidual  superelements.   Therefore,  the  residual  super- 
element  may  contain  scalar  degrees  of  freedom  representing 
linear  model  formulations, 
b.   TSTEPNL  Data 

The  TSTEPNL  data  card  description  is  provided  in 
Section  2  of  the  MSC/NASTRAN  User's  Manual.   The  input  fields 
specify  the  time  step  size,  the  number  of  steps,  and  the  out- 
put interval  as  well  as  the  nonlinear  iteration  options. 

The  nonlinear  iteration  data  are  provided  for 
controlling  the  accuracy  and  stability  of  the  nonlinear 
calculations  and  they  will  not  effect  linear  solutions.   The 
controls  are  similar  to  those  on  the  NLPARM  bulk  data  used  in 
static  nonlinear  analysis. 

As  discussed  above  two  types  of  errors  or  insta- 
bility are  encountered  in  the  transient  solution.   The  load 
iterations  are  similar  to  static  load  analysis  and  may  diverge 
at  one  or  more  time  steps  when  large  nonlinear  stiffness 
changes  occur.   Load  iteration  divergence  may  be  cured  by 
smaller  time  steps  or  by  stiffness  matrix  updates.   The 
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diagnostic  output  value  LAMBDA-I  (from  DIAG  50)  indicates 
the  rate  of  growth  of  this  error. 

The  second  type  of  divergence  occurs  in  the  time 
domain  when  the  load  iterations  do  not  converge,  or  are  limi- 
ted by  the  user.   The  errors  may  increase  for  each  new  time 
step  until  they  cause  the  code  to  abort.   This  problem  may  be 
cured  by  allowing  more  load  iterations  or  by  using  smaller 
time  steps.   The  diagnostic  output  value  LAMBDA-T  (from  DIAG 
50)  indicates  the  rate  of  growth  of  this  error. 

The  size  of  the  time  step  is  of  primary  importance 
in  transient  analysis  to  control  accuracy  and  avoid  diverging 
solutions.   The  use  of  several  subcases  with  different  time 
step  sizes  to  correspond  with  nonlinear  activity  is  recommended 
as  a  normal  procedure.   When  the  user  wishes  to  avoid  overly- 
small  time  steps  for  highly  nonlinear  zones,  he  may  increase 
the  iteration  limits  and/or  the  error  tolerance  and  the  diver- 
gence limits. 

c.   Restarts  and  Data  Base 

Solutions  66  and  99  share  the  same  data  base/ 
superelement  file  storage  logic  for  the  nonlinear  tables  and 
matrices.   Therefore,  the  restart  system  for  transient  analysis 
is  allowed  to  use  either  a  previous  static  or  transient  non- 
linear analysis  as  initial  conditions. 

For  restarting  from  previous  static  analysis  only 
the  first  subcase  is  affected.   Simply  provide  a  data  base 
created  in  SOL  66  and  specify  the  parameter. 
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PARAM,SLOPPID,N 
where   N   was  the  printed  value  of  LOOP  for  the  desired 
static  solution.   Constraint  sets,  direct  input  matrices, 
mass,  and  damping  may  be  changed. 

Restarts  from  a  previous  transient  execution  are 
available  for  a  number  of  cases.   If  the  same  model  is  to  be 
reexecuted,  only  the  residual  superelement  needs  to  be  re- 
assembled (SEMA,  SELA=0) .   If  the  end  results  for  the  previous 
transient  run  are  to  be  used  for  the  initial  conditions  at 
t=0 ,  add  the  parameter  LOOPID=N  where   N   corresponds  to  the 
Nth  subcase  of  the  previous  run  and   N   dummy  SUBCASE  commands 
to  start  the  residual  case  control  execution. 

The  normal  restart  for  a  transient  run  is  to 
continue  from  the  last  step  of  a  previous  subcase  with 
different  loads  and/or  TSTEPNL  data.   Provide  the  following 
parameters : 

LOOPID  =  N  -  Start  from  the  Nth  subcase 

Stime  =  t  -  Start  from  time  t. 

Note  that  constraint  sets  should  not  be  changed 
to  avoid  incompatible  matrix  sizes. 

D.   NUMERICAL  RESULTS 

The  progressive  deformations  of  the  plate  are  plotted 
along  the  x-axis  at  five  different  time  locations  as  shown 
in  Fig.  3.6(b).   At  the  early  time  (1.5E-4  sees),  the  response 
of  the  plate  is  in  the  elastic  region.   During  the  deformation, 
plastic  hinge  lines  are  formed  along  the  clamped  edges  and  the 
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diagonals  of  the  plate  as  indicated  in  Fig.  3.6(a).   This 
response  phenomenon  agrees  well  with  Hudson's  studies 
[Ref.  4].   As  described  in  Section  A-4  of  Chapter  3,  when 
the  plate  deflection  is  significant,  the  plate  experiences 
in-plane  stretching. 

Time  history  displacement  responses  of  Grid  points 
20101,  20201,  and  20202  are  plotted  in  Fig.  3.7.   As  shown 
in  the  figure,  for  the  short  duration  impulse,  the  maximum 
displacement  occurs  at  the  center  of  the  plate  after  termi- 
nation of  shock  loading.   The  rate  of  increase  in  displacement 
is  the  maximum  at  the  termination  of  shock  loading'  as  shown 
in  Fig.  3.8.   During  the  application  of  the  shock  loading, 
the  acceleration  level  is  high,  implying  the  high  inertia 
force  in  the  system.   At  the  termination  of  the  inertia 
force  in  the  system,  the  acceleration  stays  in  the  constant 
level  until  the  plate  displacement  reaches  the  maximum. 

Time  hisotry  Von-Mises  stress  of  QUAD4  element  213  (plate 
center)  is  plotted  in  Fig.  3.7.   The  stress  variation  is 
highly  non-linear  due  to  the  material  non-linearity. 
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Figure  3.6(a).   Plastic  Hinge  Line  Formation 
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Figure  3. 6 (b) 


Progressive  Deformation  of  the  Plate 
Along  Section  A-A. 
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IV.   CONCLUSIONS  AND  RECOMMENDATIONS 

The  state-of-the-art  literature  review  revealed  that  the 
skirt  plate  is  very  effective  to  absorb  initial  explosion 
energy  before  the  shock  loading  reaches  the  main  armor.   The 
material  of  the  armor  must  possess  the  high  hardness  to  prevent 
the  severe  damage  and  the  high  ductility  to  minimize  the 
spalling.   The  composite  materials  may  serve  more  effectively 
than  the  conventional  material  to  minimize  the  damages  against 
the  powerful  modern  projectiles. 

The  non-linear  finite  element  analysis  of  the  plate  sub- 
jected to  the  air  shock  loading  verifies  the  plastic  hinge 
line  formation  and  in-plane  stretching  inducing  snap-through 
collapse.   For  the  short  duration  shock  loading,  the  maximum 
deformation  occurs  after  the  termination  of  loading,  and  the 
velocity  reaches  the  maximum  at  the  termination  of  shock. 
The  Von-Mises  stress  is  highly  non-linear  due  to  the  material 
non-linearity.   For  the  non-linear  analysis  using  MSC/NASTRAN, 
TSTEPNL  bulk  data  card  is  found  to  be  the  most  important  card 
which  controls  the  accuracy  and  the  stability  of  non-linear 
calculation. 

To  insure  the  convergence  of  solutions,  convergence  limit 
in  the  TSTEPNL  card  in  the  bulk  data  deck  should  be  considered, 
If  a  program  was  run  with  big  convergence  limit,  the  solution 
will  be  gotten  rapidly  with  small  iteration  steps  but  it  is 
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not  accurate  and  large  time  steps  can't  be  taken.   On  the 
contrary,  if  the  program  was  run  with  small  convergence 
limit,  the  solution  will  be  accurate  but  it  is  gotten  slowly 
with  large  iteration  steps.   Thus,  it  is  expensive  and  anyway 
large  time  steps  can't  be  taken  too.   Thus,  compromise  should 
be  taken  to  run  a  large  time  steps.   For  example,  if  one  non- 
linear transient  problem  is  needed  total  700  time  steps  to  get 
the  solutions,  the  following  procedure  is  recommended.   In  the 
first  subcase,  the  program  is  iterated  with  smaller  conver- 
gence limit,  10*E-7  to  get  more  accurate  solution  and  in  the 
second  subcase,  the  program  is  run  with  limit,  10*E-5,  and  in 
the  third  subcase,  the  program  is  iterated  with  convergence 
limit,  10*E-4. 

Consequently,  the  total  700  time  step  will  be  calculated 
successively. 

The  sensitivity  studies  on  shock  period  and  amplitude  is 
suggested  to  understand  their  effects  on  plate  responses. 
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APPENDIX  A 

THE  COMMAND  TO  OPERATE  DATA  BASE 

****THE  COMMAND  TO  ERASE  THE  DATA  BASE  NAMED  "MSS . S2 770P .NODE " 

//ERASE  JOB  (2770,0252) , NASTRAN, CLASS =A 

//*MAIN  0RG=NPGVM1.2  770P 

//  EXEC  PGM=IEFER14 

//SYSPRINT  DD  SYSOUT=A 

//DD1  DD  DSN=MSS.S2770. NODE, DISP=(OLD, DELETE) 

// 

****THE  COMMAND  TO  FLOT  FROM  THE  DATABASE  NAMED  "MSS . S2770 . PPLOT" 

//LEEPLCT      JOB     ( 2770 , 0252) ,NASTRAN, CLASS =C 

//*MAIN  ORG=NPGVM1.2  770P 

//♦FORMAT  PR, DDNAME=SYSVECTR, DEST=LOCAL 

//  EXEC  NASTPLOT,PLTDSN='MSS.S2770. PPLOT' 
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APPENDIX    B 
DATA    DECK    FOR   PROGRAM 


GENERAL  DESCRIPTION  OF  OATA  DECK 

The  Input  deck  begins  with  the  required  resident  operating  system  control  cards.  The  type 
and  number  of  these  cards  win  vary  with  the  Installation.  Instructions  for  the  preparation  of 
these  control   cards  tn  given  \n  Section  7.6  of  the  MSC/NASTRAN  Application  Manual. 

The  operating  system  control  cards  are  followed  by  the  MSC/NASTRAN  Data  Oeck ,  which  consists 
of  the  following  three  sections: 

1.  Executive  Control  Oeck 

2.  Case  Control  Oeck 

3.  Bulk  Data  Oeck 

The  Executive  Control  Oeck  and  the  Case  Control  Oeck  both  have  free- field  formats.  Only 
columns  1  through  72  are  used  for  data.  Any  Information  In  columns  73  through  30  "nay  appear  In 
the  printed  echo,  but  the  data  will  not  be  used  by  the  program.  As  explained  In  Section  2.4.1, 
limited  use  Is  made  of  the  data  in  columns  73  through  30  for  the  8u1k  Oata  Oeck.  If  the  last 
character  on  a  card  1s  a  comma  (not  necessarily  In  column  72),  the  next  card  Is  a  continuation  of 
this  physical  card.  Any  number  of  continuation  cards  may  be  specified,  and  together  they  form  a 
logical   card. 

The  NASTRAN  card  1s  used  to  change  the  default  values  for  certain  operational  parameters, 
such  as  buffer  size  or  the  number  of  data  lines  printed  per  page.  More  than  one  NASTRAN  card  may 
be  present.  NASTRAN  cards  are  optional,  but.  If  present,  they  must  be  the  first  cards  of  the  data 
deck.  The  NASTRAM  card  1s  a  free-field  card  (similar  to  cards  1n  the  Executive  Control  Oeck). 
Its  format  Is  as  follows: 

NASTRAN  keyword^  •  value,   keyword2  •  value,    ... 

An  alternate  format  Is: 


NASTRAN  keywordl  ■  value 
NASTRAN  keyword2  ■  value 
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Input  Data  Card   EGRID   Grid  Point  Position 


Description:  Defines  the  position  of  a  geometric  grid  point  of  the  structural  model 


Format  and  Examde: 


1 

2 

J 

4 

5 

6 

7 

8 

9 

10 

EGRID 

10 

CP 

XI 

X2 

1 

X3 

"1 

EGRID 

2 

3 

1.0 

2.0 

\ 

3.0 

Field 


Contents 


ID  Grid  point  identification  number  (Integer  >  0) 

CP  Identification  number  of  coordinate  system  In  which  the  position  of  the  grid  point 

is  defined  (Integer  >_ 0  or  blank) 

X1.X2.X3      Position  of  the  grid  point  in  coordinate  system  CP  (Real  or  Blank) 

Remarks:   1.   The  meaning  of  XI,  XZ,  and  X3  depends  on  the  type  of  coordinate  system,  CP,  as 
follows:   (see  C3RD1  card  descriptions) 


2. 


Type 

XI 

X2 

X3 

Rectangular 
Cyl indrical 
Spherical 

X 

R 
R 

Y 

9(degrees) 
9  (degrees) 

Z 

7 

♦(degrees) 

The  EGRID  card  1s  used  1n  MSGMESH  to  define  curved  edges  and  for  other  special 
purposes  (see  Part  I,  Section  3.0).  It  1s  similar  to  a  GRID  card,  except  that  it 
will  not  be  processed  by  MSC/MASTRAN. 


3.   Entries  1n  fields  7,  8  and  9  are  not  checked  by  MSGMESH. 
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Input  Data  Card   GRIDG   Linear  and  Higher  Order  Interpolation 


Description :  Defines  a  field  of  grid  points  and  generates  GRID  cards 


Format  and  Example: 


1 

2 

3 

4 

5 

5 

7 

8 

9 

10 

GRIOG 

FID 

CD/ CI 

PS 

I 

GA 

G8 

GC 

SEID 

GRIDG 

10 

+0 

18 

-1 

-2 

-3 

L 

M 

m — 

GD 

N 

GE 

GF 

GG 

' 

GH 

: 

CI 

♦€2 

o 

-4 

1 
|*G2 

EKAB) 

E2(A8) 

El(BC) 

E2(5C) 

EKCD) 

E2(CD) 

El  (DA) 

E2(DA) 

+G2 

1 

6 

5 

8 

7 

EKEF) 

E2(EF) 

El(FG) 

E2(FG) 

El(GH) 

E2(GH) 

EKHE) 

E2(HE) 

1 

El(AE) 

E2(AE) 

E113F) 

E2(BF) 

EKC3) 

E2(CG) 

EllDH) 

E2(DH) 

| 

Field  Contents 

FID         ID  number  of  a  GRIDG  field  (Integer,  1  _<_  FID  £99).   Must  be  unique  for  all  GRIDG 
cards.  When  entry  Is  negative.  Internal  EGRID  cards  are  generated. 

CD/CI        Prescribes   displacement   and/or   interpolation  coordinate   systems   (Integer  or 
blank).  See  Remarks  4  and  5. 

PS         PS  entry  on  generated  GRID  cards  (any  of  the  digits  1-6  with  no  Imbedded  blanks) 
(Integer  >_  0  or  blank) 

L,M,N        Grid  field  dimensions  when  an  entry,  L,  M,  and/or  H    is  positive.   (See  Figure  3.) 
Absolute  value  is  a  LID  entry  on  a  LIST  Input  card,  when  an  entry  is  negative. 

SEID        SEID  entry  on  generated  GRID  cards  ( Integer  >_0   or  blank) 

GA,...,GH    Grid  point  Identification  numbers,  GIDs,  of  grid  field  vertex  points  (Integer  /  0  or 
blank.).  See  Remarks  3  and  5. 

CI         CP  entry  on  generated  GRID  cards'  (Integer  >_  0  or  blank),  and  grid  field  inter- 
polation coordinate  system  definition.  See  Remarks  4  and  5. 

El(XY)       ID  number  of  edge  grid  point  (closer  to  X)  on  side  XY  (Integer  A  0  or  blank).  See 
Remarks  3  and  6. 

E2(XY)       ID  number  of  edge  grid  point  (closer  to  Y)  on  side  XY  (Integer  r  0  or  blank).  See 
Remarks  3  and  5. 

(Continued) 
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GRIOG  (Cont.) 

Remarks:   1.    See  Part  I,  Sections  2  and  3  for  discussions  on  the  GRIOG  Input  card. 
2.   Required  data  by  type  of  field 


Type 

Required  Data 

LINE 

TRIA 
QUAD 
HEX 

FID,  L,  GA,  GB 
FID,  L,  GA,  G3,  GC 
FID,  L.  H,  GA,  GB,  GC ,  GD 
FID,  I,  M,  N,  GA  thru  GH 

Only   the   parent 
Interpolation. 


and   first   continuation   card   are   applicable   for   linear 


Vertex  and  edge  points  an  defined  on  GRID  or  EGRID  input  cards,  or  on  GRID  cards 
generated  by  previously  processed  GRIDG  Input  cards.  Field  vertex  points  are 
automatically  equivalenced  to  corresponding  entries  GA,  ...,  GH,  when  these  entries 
are  positive.  GRID  cards  are  generated  only  for  active  field  grid  points.  See 
Part  I,  Section  10. 1. 


CO /CI 


the  entry   is 

The  CI  entry 

exist  and  an 


entry  may  be  used  to  prescribe  both  the  displacement  coordinate  system,  CO, 
and  the  Interpolation  coordinate  system,  CI.   3oth  CO  and  CI  ar^    prescribed  when 
the  entry  1s  positively  signed.    Only  CD  Is  prescribed  when 
unsigned.   Only  CI  is  prescribed  when  the  entry  is  negatively  signed. 
on  the  first  continuation  card  may  be  used  when  both  CO  and  CI 
different. 

CI  references  a  CID  entry  on  a  coordinate  system  definition  Input  card.  When  CI  is 
entered  as  0,  interpolation  is  in  basic  rectangular  coordinates.  When  CI  is  not 
entered,  interpolation  1s  in  the  coordinate  system  that  positions  the  field  vertex 
and  edge  points.  In  this  case,  all  vertex  and  edge  points  must  be  defined  in  the 
same  coordinate  system  (CP  entries  on  GRID  and/or  EGRID  cards).  See  Part 
Section  3.1  and  Table  2. 


I, 


The  $  interpolation   point 

the  9  interpolation  point  components, 

offset  by  a  t360°  by  appending 


components,   in   spherical   coordinate   systems,   or 
ents,  1n  cylindrical  coordinate  systems,  may 
a  tS  or  a  tC  to  grid  point  entries. 


be 


7.  If  in  entry  for  GO  is  present  in  a  TRIA  field,  grid  point  D  replaces  grid  point  C 
in  Interpolating  grid  point  positions  from  side  CA.  Grid  point  0  is  usually 
defined  by  in  EGRID  card.     M  must  be  blank.      See  Part    I,    Section  3.5. 

3.  If  both  El(XY)  and  E2(XY)  are  entered,  the  positions  of  generated  grid  points  an 
determined  by  cubic  Interpolation  from  edge  XY.  If  E2(XY)  is  blank,  positions  an 
determined  by  quadratic  Interpolation.  If  both  El(XY)  and  E2(XY)  are  blank, 
positions  an  determined  by  linear  interpolation.  Either  none  or  one  edge  point, 
E(XY)  may  be  specified  for  each  edge  of  a  TRIA  field.  El  (CO )  positions  a  grid 
point  on  edge  CA  of   TRIA  fields.      See  Part    I,    Section   3.7. 

9.        Intermediate  continuation   cards   cannot   be  omitted  even    if  all    entries   are   blank.- 
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Input  Data  Card   CGEN   Generate  Element  Connection  Bulk  Data  Cards,  Form  1 
Description:  Generate  element  connection  cards  for  element  types  listed  under  Remark  1 


Format  and  Example: 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

CGEN 

TYPE 

FEID 

PID 

PID 

DIR 

TH 

EIDL 

EIDH 

C3EN 

QUA04 

100 

10 

1 

+C2 

IT 


♦<2 


TA 


TC 


0.002 


0.002 


0.001 


TD 


0.001 


Field 
TYPE 

FEID 

PID 

FID 
DIR 

TH 


EIDL 

EIDH 
TA,T8,TC,TD 


Contents 

Any  of  the  mnemonics  listed  under  Remark  1 

The  absolute  value  of  FEID  is  the  element  Identification  number,  EID,  of  first 
element  1n  element  set  (Integer  **  0).     See   Remark  4. 

ID  of  an  element  property  card  or  PGEN  card  (Integer  >  0  or  blank).  Must  be  unique 
with  respect  to  all    other  property  card   10s  and  PGEN   IDs.      See   Remark   5. 

ID  of  associated  grid  point   field   defined  on  a  GRIDG  or  GRIDU  card 

One  of  the  entries  L,  M,  or  V  (Default  entry  1s  L).  Direction  of  elements  In  the 
field.     See  Table  4  of  Part   I,   Section  7.1. 

Material  prooerty  orientation  specification  (Real  or  Blank;  or  Integer  >_ 0  and  less 
than  1,000,000).  If  Real  or  31ank,  specifies  the  material  property  orientation 
angle  In  degrees.  If  Integer,  the  orientation  of  the  material  x-axis  is  along  the 
projection  onto  the  plane  of  the  element  of  the  x-axis  of  the  co-ordinate  system 
specified  by  the  Integer  value. 

ID  of  element,  below  which  no  element  connection  cards  are  generated  (Integer  >  0 
or  blank) 

ID  of  element,  above  which  no  element  connection  cards  ire  generated  (Integer  >  0 
or  blank) 

Thickness  of  corner  points  of  QUAD  and  TRIA  fields  when  TYPE  is  TRIA3,  TRIA6,  0UA04 
or  QUAD8   (Real   or  Blank).     See  Remark  6. 


Remarks:     1, 


Element  types:  R0O,  PL3TEL,  CFTUBE,  TUBE,  TRIA3,  TRAPRG,  TRIARG,  TRIA6,  TRIM6, 
TRIAX6,  QUAD4,  SHEAR,  QUAD8,  HEXA1,  HEXA2,  HEXA8,  HE*A20  and  HEX20.  The  mnemonics 
HEXA8   and  HEXA20  refer  to  HEXA  elements  with  8  or  9   through  20   nodes,    respectively. 


(Continued) 
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CGEM   (Cont.) 


2.  Dements  are  tabulated  by  class  In  Table  3,  of  Part  I,  Section  7.0.  Line  class 
elements  may  be  generated  1n  all  field  types.  Tr1a3  and  Tr1a6  class  elements  may 
be  generated  1n  TRIA  and  0UA0  fields.  Ouadl  and  Quad9  class  elements  may  be 
generated  1n  TRIA,  QUAD  and  HEX  fields.  Only  QUAD*  and  QUAC8  element  tyoes  are 
valid  1n  TRIA  fields.  TRIA3  and  TRIA6  elements  are  generated,  1n  this  case,  along 
the  field  edge  corresponding  to  the  element  direction.  Hex8  and  Hex20  class 
elements  may  only   be  generated   1n  HEX   fields. 

3.  MSGMESH  automatically  generates  PENTA  elements  in  place  of  HEXA  and  HEX20  elements, 
wnen  adjacent  element  edges  would  have  the  same  grid  point  Identifiers.  TRIA6 
elements  replace  QUAD8  elements  and  TRIA3  elements  replace  0UAC4  elements,  when 
adjacent  element  vertices  would  have  the  same  grid  point  Identifiers.  Otherwise, 
MSGMESH  suooresses  the  generation  of  an  element  which  would  reference  Identical 
grid  point  Identifiers.  Elements  are  generated  only  when  all  vertex  points 
exist.     Deleted  mldslde  nodes   suopress  the  generation  fo  TRIMS  and  TRIAX6  elements. 

4.  EIDs  are  Incremented  when  the  FEID  entry  1s  positive  and  decremented  when  the  FEID 
entry   1s  negative. 

5.  The  P ID  entry  may  be  blank  for  those  elements  whose  element  connection  card  permit 
a  blank  property   ID. 

6.  Element  thicknesses  may  be  generated  for  TRIA3,  TRIA6,  CUA04  and  CUA08  elements 
located  In  TRIA  and  QUAD  fields.  Interpolation  from  the  thicknesses  at  the  field 
vertices  Is  based  on  field  topologies  with  equal  spadngs  of  thicknesses  for  both 
GRIOG  and  GRIDJJ  fields.  Unequal  spadngs  of  thicknesses  are  defined  for  GRIDG 
fields   specified  with  unequal    spadngs  of  grid  points. 
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Inout   Oata   Card       MAT!        Materia)    Property  Definition,    Form   1 

Description:        Defines     the    material     properties     for     linear,     temperature-independent,     isotropic 

materials 


Format    and   Example: 


1 

2 

3 

4 

5 

6 

7 

3 

9 

10 

MAT1 

MID 

r 

G 

Ml! 

9H0 

A 

TREF 

GE 

| 

MAT! 

17 

3. +7 

0.33 

4.28 

6.5-6 

5.37*2 

0.23 

ABC 

i           ST 

SC                    SS 

MCSID 

1 

*8C 

20. *4 

IS.*4 

12.+4 

1003 

Field 

MID 

E 

G 

NU 

RH0 

A 

TREF 

GE 

ST.SC.SS 

MCSID 


Contents 

Material    identification  number   (Integer  >   0) 

Young's  modulus    (Real    or  blank) 

Shear  modulus    (Real    or  blanic) 

Polsson's   ratio    (-1.0   <  Real  _<  0.5  or  blank) 

Mass   density    (Real) 

Thermal    expansion   coefficient    (Real) 

Thermal    expansion   reference  temoerature    (Real) 

Structural    element   damping  coefficient    (Real) 

Stress  limits  for  tension,  compression,  and  shear  (Real).  (Used  only  to  compute 
margins  of  safety  in  certain  elements;  they  have  no  effect  on  the  computational 
procedures. ) 

Material   Coordinate  System  identification  number   (Integer   >  0  or  blank) 


Remarks:      1.       The   material    identification   number   must  be  unique   for  all    MAT1,   MAT2,   MAT3  and  MAT9 
cards. 

2.  MAT1  materials  may   be  made   temperature  dependent   by   use  of  the  MATT1  card. 

3.  The     mass     density,     RH0,     will     be     used     to     automatically     compute     mass     for     all 
structural   elements. 

4.  Weight   density   may   be  used  in  field  5  if  the  value   l/g   is   entered  on   the  param  card 
WTMASS,  where  g  is  the  acceleration  of  gravity    (see  Section  3.1.5). 

5.  MCSID   must    be    nonzero    if    the   CURV   module    is    used    to   calculate    stresses    or   strains 
at   grid  points  on  plate  and  shell   elements  only. 

6.  To  obtain   the   damolng   coefficient,   GE,   multiply  the  critical    damping   ratio  C/Ca,    by 
2.0. 

(Continued) 
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MAT1  (Cont.) 


7.  Either  E  or  S  must  be  specified  (I.e.,  nonblank). 

8.  If  any  one  of  E.  G,  op  NU  is  blank,  it  *i  1 1  be  computed  to  satisfy  the  identity 
E  ■  2(1+NU)G;  otherwise,  values  supplied  by  the  user  will  be  used.  This 
calculation  is  only  made  for  initial  values  of  E,  3,  and  NU . 

9.  If  E  and  NU  or  G  and  NU  are  both  blank,  they  will  both  be  given  the  /alue  0.0. 

10.  Implausible  data  on  one  or  more  MAT1  cards  will  result  in  a  warning  message. 
Implausible  data   1s   defined   as   any   of   E  <  0.0,  G  <  0.0,  0.5  <  NU  <  0.0,   or 


1   - 


Z( 1*NU^G    I    *   0*°*    (exceQt    for  case*   covered   by   Remark    9). 
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Input  Data  Card   MATS!   Material  Stress  Dependence 

Description:  Specifies  table  references  and  material  properties  which  are  stress-dependent  for 
use  m  -naterlal  nonl  1neaM ty  applications.  This  card  will  be  activated  If  MAT1  with  the  same  MID 
1s  being  used  1n  the  material  nonlinear  solution  sequence  (66). 


Format  and  Example: 


10 


MATS1 

MID 

Rl 

TYPE 

3 

YF 

HR 

LIMIT1    I^XCj 

MATS1 

17 

28 

PLASTIC 

2 

3 

2.*4 

+ABC 

LIMITZ 

+ABC 

Field 

MID 
Rl 

TfPt 

3 


YF 


HR 


UMIT1 
LIMITZ 


Contents 

Identification  number  of  a  MAT1  card  (Integer  >  0) 

Identification  number  of  a  TABLES1  card  (Integer  >  0) 

Type  of  material  nonllnearlty  (BCD)  MLELAST  (Nonlinear  elastic)  or  PLASTIC 
(Plasticity  theory) 

Wort  hardening  slope  (slope  of  stress  vs.  plastic  strain)  1n  units  of  stress 
(Real).  For  elastlc-perfectly  plastic  cases,  B  a  0.0.  For  more  than  a  single 
slope  in  the  plastic  range,  the  stress-strain  data  must  be  supplied  1n  a  TABLES1 
card  referenced  In  field  3.     See  Remarks   1  and  2. 

Yield  function  (Integer);  one  of  the  following  values:      (Default  »  1) 

1  (von  M1ses) 

2  (Tresca) 

3  (Mohr-Coulomb) 

4  (Orucker-Prager) 

Hardening  Rule  (Integer);  one  of  the  following  values:  (Default  ■  I) 

1  (Isotropic) 

2  (Kinematic) 

3  (Combined  Isotropic  and  cinematic  hardening) 

Parameter  used  in   the  yield  function  specification  (Real).  See  Remark  4. 

Parameter  used  in  the  yield  function  specification  for  yield  functions  3  and  4 
(Real ).  See  Remark  4. 


33 


MATS1  (Cont.) 


Remarks:  1.  If  type  •  NLELAST,  the  stress-strain  data  given  fn  the  TA8LES1  card  will  be  use- 
determine  the  stress  for  a  given  value  of  strain.  The  values  8,  YF,  HR,  LIM 
and  LIMIT2  will  not  be  used  1n  this  case. 

2.  If  type  »  PLASTIC,  either  the  table  Identification,  Rl,  or  the  wort  harder 
slope,  3,  may  be  given  but  not  both.  If  the  table  ID  Is  omitted,  the  \ 
hardening  slope,  B,  specified  in  M-eia  3  1s  defined  as 


f 


r 


or 


V 


1  * 


where   E    1s    the   elastic   modulus   and   f*    1s   the   slope  of   the  uniaxial    stress-str 
curve  1rr  the  plastic  region. 

If  Rl   Is  given,   TABLES!  entries  of  stress-strain  data  must  conform  to  the  follow 
rules: 

a.  If  type  ■  NLELAST,   the  curve  must  pass  through  origin  (x.,*0.,  yfa0.,   for  s 
value  of  1).     This  option  1s  valid  only  for  line  elements. 

b.  If  type  »  PLASTIC,  only  data  for  the  first  quadrant  need  be  supplied.  Th 
are  two  options  available  to  specify  the  data.  In  case  L  (defined  as  F9RI 
In  the  TABLES1  card)  the  first  point  1s  at  the  origin,  the  second  point  Is 
the  yield  stress  given  1n  field  3  of  the  MATS1  card.  The  slope  of  the  1 
joining  the  origin  to  the  yield  stress  must  be  eaual  to  the  value  of  E  g1' 
on  the  MAT1  3u1  k  Oata  card.  In  case  2  (defined  as  FORM-1  in  the  TA8L1 
card)  the  slooes  of  the  curve  1n  the  plastic  range  in  supplied.  See  I 
TA8LE51  card  description  for  additional   details. 

LIMIT1  and  LIMIT2  are  parameters  used  1n  the  yield  function  definition  as  follows 


Yield  Function 

LIMIT1 

LIMIT2 

von  M1ses  (1) 
or  Tresca  (2) 

Mohr-Coulomb  (3) 

or  Orucker-Prager  (4) 

Yield  stress  1n  tension,  a 

y 

Cohesion, c  [in   stress  units) 

Mot  used 

Angle  of  internal 
friction  s  {\n   degree; 

See     Section     2.14.1.2     \n     the    MSC/NASTRAN     Application     Manual     for    a     detail 
description  of  yield  function  and  hardening  rules. 
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MATS1  (Cont.) 


*-   £ 


NLELAST 


0. 


PLASTIC 
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Input  Data  Card       ~ABL£S1       Material    Property   Table,   Form   1 

Description:       Defines    a    tabular    function    for    stress-dependent    material     properties    sucn    as    the 
stress-strain   curve  and   creep   parameters 


Format    and   Examole: 


1 

2 

3 

4 

s 

6 

7 

a 

9 

10 

:abl£si 

10 

I 

X 

^><3]><C1^><CJ 

TABLESl 

32 

I 

A8C 

xl 

*l 

x2 

y2 

x3 

*3 

*4 

,  — 
*4 

*8C 

0.0 

0.0 

.01 

10000. 

.02 

15000. 

EN0T 

(etc.) 


Field 

10 

I 


Contents 
Table  Identification  number  (Integer  >  0). 
Not  used 
Tabular  entries  (Real) 


Remarks :   I.   The  x^  must  be  1n  either  ascending  or  descending  order  but  not  both. 

2.  Jumps  between  two  points  (x^  *  xi*i)  are  allowed,  but  not  at  the  end  points. 

3.  At  least  two  entries  must  be  present. 

4.  Any  x-y  entry  may  be  ignored  by  placing  the  BCD  string  SKIP  in  either  of  the  two 
fields  used  for  that  entry. 

5.  The  end  of  the  table  is  indicated  by  the  existence  of  the  3C0  string  EN0T  in  eitner 
of  the  two  fields  following  the  last  entry.  An  error  is  detected  if  any 
continuation  cards  follow  the  card  containing  the  end-of-table  Hag  EN0T. 

5.   TABLESl  is  used  to  input  a  curve  in  the  form  of 


iM 


*nere  X  is  input  to  the  table  and  Y  is  returned.  The  table  look-uo  yj(x),  x  »  X, 
1s  performed  using  linear  interpolation  within  the  table  and  linear  extrapolation 
outside  the  table  using  the  last  two  end  points  at  the  appropriate  table  end.  At 
Jump  points  the  average  y-r(x)  is  used. 
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Inout  Data  Card       P13AQ2       Pressure  Load  on  a  Two-Dimensional    Structural   Element 

Descriotion:      Defines    a    uniform   static    pressure    load    applied   to    two-dimensional    elements.      Only 
QUAD4,   5HE.AR,   TRIA3,   or  TRIAX6  elements  may  have  a  pressure  load  applied  to   them  via   this  card. 

Format  and  Example: 


1 

2 

3 

4 

5 

6 

7 

a 

9 

10 

PL3AD2 

SID 

0 

EIO 

EIO 

EID 

EID 

EID 

EID 

PL3AD2 

21 

-3.5 

4 

16 

1    2 

|         i 

Alternate  Form: 


PL0AO2 


SID 


EID1 


"THRU" 


EID2 


PL3AD2 


30.4 


15 


THRU 


48 


Field 
SID    . 


EIO  ) 
EID1> 
EI02) 


Contents 
Load  set  identification  numoer  (Integer  >  0) 
Pressure  value  (Real ) 

Element  Identification  number  (Integer  >  0;   EID1   <   EID2) 


Remarks:      1.       EIO  must  be  0  or  blank   for  omitted  entries. 

2.  Load    sets    must    be    selected    in    the    Case    Control    Decic    (L3AD=SID)    to    be    used    by 
MSC/NASTRAN. 

3.  At  least  one  positive  EID  must  be  present  on  each  PL3A02  card. 

4.  If     the     alternate     form    is    used,     all     elements     EID1     through     EID2    must    be     two- 
dimensional  . 

5.  The    direction   of    the   pressure    is    computed   according    to    the    right-*and   rule   using 
the  grid  point  sequence  specified  on  the  element  card.     Refer  to  the  PL3AD  card. 

5.   All  elements  referenced  must  exist. 

7.   Continuation  cards  are  not  allowed. 
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[nput  Data  Card   NLPARH   Parameters  for  Nonlinear  Analysis  Control 

Pescri  ption:   Defines  a  set  of  parameters  for  nonlinear  analysis  iteration  strategy 


Format  and  Example: 


10 


NIP ARM 

ID 

INC 

or 

<METH0D 

<S7EP 

MAX  ITER 

C3NV 

INT3UT 

NLPARM 

15 

5 

AUT30N 

PW 

-NL1 

EPSU 

EPSP 

EPSW 

MAXOIV 

MAXqN 

MAXLS 

FSTRES3 

*NL1 

1.0E-3 

1.0-4 

1.0-5 

1 

20 

3 

.02 

Field 

ID 

INC 

OT 

KMETH0P 

KSTEP 

V63   MAX  ITER 
C0NV 

INT3UT 

EPSU, EPSP, 
EPSW 

MAXOIV 

MAXQN 

MAXLS 

FSTRESS 


Contents 

Identification  Number  (Integer  >  0) 

Number  of  increments  (Integer  >  0) 

Incremental  time  Interval  for  creep  analysis  (Real  _>.  0)  (Default  *  0.;  no  creep) 

Method  for  controlling  tangent  stiffness  updates  (3CD  -  AUT0, ITER, SEMI ,LSON,AUT0CN, 
or  SEMIQN)  (No  default).  See  Remark  5. 

Number  of  iterations  before  a  .matrix  update.    Used  only  for  ITER  and  LSQN 
methods.   ( Integer  >_  0  or  blank)  (Default  »  1) 

Limit  on  total  iterations  for  each  load  increment  (Integer  k   0)  (Default  ■  20) 

Flags  to  select  convergence  tests  (9CD  ■  U,  P.  W,  or  any  combination  of  the 
letters)  (Oefault  *  W).   See  Remarks  3  and  4. 

Intermediate  output  flag  (BCD  »  YES  or  N0)  (Oefault  »  YES).   See  Remark  5. 

Error  tolerance  limits  for  determining  convergence  (Real  >  0.0)  (Defaults  ■  l.DE-2, 
1.0E-3,  1.0E-4,  respectively) 

Limit  on  diverging  Iterations  (Integer  >  0)  (Oefault  *  1) 

Maximum  number  of  Quasi -Newton  correction  vectors  to  be  saved  on  the  data  base 
(Integer  _>  0)  (Default  =»  10) 

Maximum  number  of  linear  search  operations  per  vector  iteration   (Integer  >  0) 
(Default  «  3) 

Fraction  of  effective  stress  (o)   used  to  determine  the  subincrement  size  in  the 
material  routines  (0.  <  Real  <  1.0)  (Default  »  0.15).   See  Remark  9. 


(Continued) 
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NLPARM  (Cont.) 


Remarks:   1.   The  NLPARM  3u!k  Data  card  is  requested  by  the  Case  Control  card  NLPARM  =  ID.   EacM 
solution  subcase  requires  a  loading  condition  and  an  NLPARM  request. 

2.  In  case  of  static  analysis  (OT  *  blank),  INC  is  the  numoer  of  equal  subdivisions  of 
the  net  load  change  defined  for  the  SUBCASE.  Aoplied  loads,  gravity  loads, 
temperature  sets,  enforced  displacements,  etc.,  define  the  new  loading  condition. 
The  differences  from  the  previous  case  are  divided  by  INC  to  define  the  incremental 
values.  In  case  of  creep  analysis  (OT  >  0.),  INC  is  the  number  of  time  step 
Increments. 

3.  The  test  flags  (U,P,'<0  and  the  error  limits  (EPSU,  etc.)  define  the  convergence 
criteria  (U  »  Displacement  error  test,  P  »  Load  equilibrium  error  test,  A  =  «or< 
error  test).  If  the  internally  calculated  error  fractions  are  less  than  the 
requested  limits,  the  iteration  stoos,  the  results  ire  processed,  and  the  orogrjm 
continues  to  the  next  load  increment  or  subcase.  All  requested  error  criteria  must 
be  satisfied  for  a  "convergent  solution". 

*.  If  convergence  does  not  occur  for  a  particular  load  increment,  the  number  of 
iterations  is  limited  to  MAXITER.  If  MAXITER  is  a  positive  number,  the  results 
will  be  processed  and  the  program  will  continue  with  the  next  load  increment.  If 
MAXITER  is  a  negative  value,  the  solution  will  terminate  if  convergence  is  not 
achieved. 

5.  The  basic  nonlinear  solution  methods  are  the  "QFGS"  Quasi-Newton  and  modified 
"Newton-Rapnson"  processes  whereby  the  out-of-baiance  nonlinear  loaas  are  neasureo 
and  recycled  to  solve  for  improved  displacements.  The  convergence  of  this 
iteration  process  is  generally  improved  by  either  performing  "line  searches"  or  oy 
frequent  updates  to  the  "tangent  stiffness  matrix"  to  account  for  current  plastic 
and  geometric  effects.  However,  the  matrix  updates  may  be  relatively  costly  jnd 
may  be  detrimental  if  the  estimates  are  far  from  the  correct  answer.  Thus,  the 
user  is  given  the  following  KMETH0O  options  to  provide  control  over  a  large  range 
of  problem  types. 

a.  If  the  AUT8  option  is  selected,  the  program  will  automatically  select  the  most 
efficient  strategy  based  upon  convergence  rates.  At  eacn  step  the  number  of 
iterations  required  to  converge  is  estimated.  If  this  exceeds  the  MAXITER 
limit,  or  if  the  CPU  time  exceeds  the  time  to  perform  a  stiffness  matrix 
update,  the  module  will  perform  a  stiffness  update.  If  divergence  occurs  for 
MAXDIV  successive  iterations,  a  stiffness  update  is  also  performed. 

b.  If  the  SEMI  option  is  selected,  for  each  load  increment  the  program  will:  (1) 
perform  a  single  displacement  iteration  based  upon  the  new  load;  (2)  update 
tne  stiffness  matrix  for  the  estimated  geometry  and  material  properties;  and 
(3)  resume  the  normal  AUTO  iteration  method. 

c.  If  the  ITER  option  is  selected,  the  stiffness  matrix  will  be  updated  only 
after  eacn  KSTE?  iteration.  The  count  is  reset  to  zero  for  each  new  matrix  or 
load  step. 

d.  If  one  of  the  AUT0QN,  SEMIQN,  or  the  LSQN  ootions  is  selected  the  "Duasi- 
Newton"  and  "Line  Search"  procedures  are  used.  3asically  the  line  search  (LS) 
will  scale  the  displacement  increment  vector  to  minimize  the  error  at  each 
iteration.  If  the  line  search  provides  a  substantial  improvement,  the  Duasi- 
Newton  (QN)  operations  use  the  local  tangent  of  the  error  to  correct  the 
tangent  matrix  solution  for  subsequent  iterations.  Otherwise  the  operations 
behave  much  like  the  AUT0,  SEMI,  and  ITER  options. 

(Continued) 
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NLPARM  (Cont.) 


6.  Output  requests  for  ELF0RCE  and  STRESS  made  In  Case  Control  will  be  processed  if 
INT0UT  -  YES  (tne  default  value).  If  INT9UT  -  H9,  these  output  requests  are 
processed  for  only  the  last  load  factor  In  a  subcase.  Requests  for  OISP  printout 
and  structure  plots  in  also  processed  during  the  nonlinear  Iteration  1f  PARAM, 
NLOISP,  1  1s  provided  In  the  Bulk  Data  deck. 

7.  If  a  diverging  solution  1s  encountered  for  MAXOIV  successful  Iterations,  the  AUT0, 
AUT3QN,  SEMI,  and  SEMIQN  options  -ill  perform  a  tangent  stiffness  update.  If, 
after  the  stiffness  update,  the  solution  Is  still  diverging,  the  program 
terminates.  For  ITER  and  LSQN  options,  diverging  cases  cause  immediate 
termination. 

8.  The  unit  of  the  Incremental  time  Interval  under  OT  field  mist  be  consistent  with 
the  unit  used  on  the  CREEP  card  to  define  the  creep  characteristics.  Total  creep 
time  for  the  subcase  may  be  found  from  OT  multiplied  by  the  value  in  the  INC. 

9.  The  number  of  sublncrements  1n  the  material  routines  (elastoplastlc  and  creep)  are 
determined  such  that  the  subincrement  size  Is  approximately  FSTRESS  *  o. 

FSTRESS  1s  also  used  to  establish  a  tolerance  for  elastoplastlc  material  in  the 
material  routine,  I.e., 


<5  1n  yield  function  <  FSTRESS  *  yield  stress 


If  the  limit  is  exceeded  at  the  converging  state,  the  program  will  exit  with  a 
fatal  error  message.  Otherwise  the  stress  state  is  adjusted  to  the  current  yield 
surface,  resulting  in  S  »  0. 
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Input  Data  Card       SPCG       Constraint  Set  Generator,    Standard  Form 

(Description:     Generates  SPC1  cards   that  constrain  a  subset  of  grid  points   1n  a  grid  point   field 


Format  and  Example: 


1 

2 

3 

4 

5 

6 

7 

a 

9 

10 

SPCG 

SID 

FID 

C 

01 

D2 

SPCG 

15 

3 

246 

0002 

0312 

Field 

SID 

FID 

C 

01,  D2 


Contents 

ID  number  of  single-point  constraint  set  (Integer  >  0) 

ID  "number  of  a  GRIDG  or  GRIDU   field 

Component   numbers    (any   combination   of   the   digits    1-6   with   no    Imbedded   blanks)    1n 
the  displacement  coordinate  system  (CD  entry   on  GRIDG  card) 

Location   number  pair  of  order- Independent  pivot  points  that  define  a  subset  of  grid 
points   In  FID.     See  Part   I,   Section  2.5. 


PemarVs:     1.       See  Part   I,   Section  5.0   for  a  discussion  on  the  SPCG  Input  card. 

2.       SPC    entries    on    generated    SPC1    cards    are    made    only    for    active    grid    points.      See 
Part   I,   Section  10.1. 
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Input  Oata  Card   DL3AQ   Dynamic  Load  Comoinatlon  (Superposition) 

Oescri  pt 1  on :  Defines  a  dynamic  loading  condition  for  frequency  response  or  transient  response 
proolems  as  a  linear  combination  of  load  sets  defined  via  RL9AD1  or  RL3A02  cards  (for  frequency 
response)  or  TL3AD1  or  TL3AD2  cards  (for  transient  response) 

Format  and  Examples: 


1 

2 

3 

4 

S 

i 

7 

1 

9 

10 

DL3AD 

SIQ 

S 

SI 

LI 

S2 

L2 

S3 

L3 

)UWO 

17 

1.0 

2.0 

6 

-2.0 

7 

2.0 

3 

-A 

S4 

L4 

-etc«- 

*A 

-2.0 

9 

-etc- 


F'eld 

SID 

S 

Si 

Li 


Contents 

Load  set  identification  number  (Integer  >  0) 

Scale  factor  (Real ) 

Scale  Factors  (Real  ) 

Load   set   identfi cation   numbers   defined   via   card   types   enumerated   above 
(Integer  >  0) 


Remarks :   1.   The  load  vector  being  defined  by  this  card  is  given  by 


W-$[  5,1^,1   • 


2. 

3. 
4. 

5. 

6. 

7. 
8. 
9. 


The  Li  must  be  unique. 

SID  must  be  unique  from  all  Li. 

Nonlinear  transient  loads  may  not  be  included;  they  ir9    selected  separately  in  the 
Case  Control  Oeck. 

Linear  load  sets  must  be  selected  in  the  Case  Control  Deck  (DL3AD=SID)  to  be  used 
by  MSC/NASTRAN. 

A  0L3A0  card  may  not  reference  a  set  identification  number  defined  by  another  DL3A0 
card. 

TL0AO1  and  TL0AO2  loads  may  be  combined  only  througn  the  use  of  the  0L3A0  card. 

RL3A01  and  R13A02  loads  may  be  combined  only  through  the  use  of  the  DL3AD  card. 

SID  must  be  unique  for  all  TL0AO1,  TL3AD2,  RL0AD1  and  RL0AD2  cards. 
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Input   Data   Card       TLJADl        Transient   Response  Dynamic  Load,   Form   1 


Description:      Defines   a   time   dependent    dynamic    load   or   enforced  motion   of    the   form 


[P(t)}    -    (A   F(t    -   t)} 


for  use  1n  a  transient  response  problem 


Format  and  Example: 


1 

2 

3 

4 

5 

6 

7 

a 

9 

10 

TL0AO1 

SID 

I 

n 

TYPE 

TID 

TLJADl 

5 

7 

9 

0 

13 

Field 
SID 

L 

M 
TYPE 

TID 


Contents 

Set  identification  number  (Integer  >  0) 

Identification  number  of  DAREA  card  set  or  a  static  load  card  set  or  a  thermal  load 
set  (In  heat  transfer  analysis)  which  defines  A  (Integer  >  0) 

Identification  number  of  DELAY  card  set  which  defines  r  (Integer  2.  D) 

Defines  the  nature  of  the  dynamic  excitation  (Integer  0,  1,  2,  3  or  blank).   (See 
Remark  2.) 

Identification  number  of  TA8LED1  card  which  gives  F(t  -  tj  (Integer  >  0) 


Remarks:   L.   If  H  1s  zero,  t  will  be  zero. 

2.   The  nature  of  the  dynamic  excitation  Is  defined  in  field  5  in  accordance  with  the 
following  table. 


Integer 

Excitation  Function 

0  or  blank 

1 

2 

3 

Force  or  Moment 
Enforced  Displacement 
Enforced  Velocity 
Enforced  Acceleration 

See  Section  2.7  of  the  MSC/NASTRAN  Application  Manual  regarding  the  use  of  the 
enforced  motion  options.  Note  that  "large  masses"  must  be  used  for  enforced 
motion.  For  heat  transfer  problems.  Field  5  must  be  blank. 

(Continued) 
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TL3A01    (Cont.) 


3.  Dynamic  load  sets  must  be  selected  1n  the  Case  Control  Oectc  {013A0=S1D)  to  be  used 
by   1SC/NASTRAN. 

4.  T10AO1  loads  may  oe  comoined  with  TL3A02  loads  only  oy  specification  on  a  DL3A0 
card.  That  Is,  the  SID  on  a  TL3A01  C3l-<j  may  not  be  the  same  as  that  on  a  TL^AOZ 
card. 

5.  SID  must   be  unique   for  all    TIJA01,   TL0AO2,   3L3A01   and  913AD2  cards. 

6.  Field  3  may  reference  sets  containing  QHBOY,  QB0Y1 ,  QBDY2,  QVECT,  and  QV0L  cards 
when  using  the  heat   transfer  option. 

7.  If  the  heat  transfer  option  is  used,  the  referenced  QVECT  data  card  may  also 
contain   references   to   functions   of   time,   and  therefore  A  may   be  a   function   of    time. 

8.  Fourier  analysis  will  be  used  1f  this  Is  selected  in  an  aeroelastic  response 
problem. 
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Inout  Data  Card       LSEQ       Static  Load  Set  Definition 

Description:      Defines   a    sequence   of   static    load   sets   to   be   applied   to    the   structural    nodel  .      The 
load   sets  may  be  referenced  by  dynamic  load  cards. 

format  and  Example: 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

LSEQ 

SID 

OAREA 

LID 

TID 

LSEQ 

100 

200 

1000 

1001 

Field 
SIO 
OAREA 
LID 

TID 


Contents 

Set  identification  of  the  set  of  LSEQ  cards     (Integer  >  0) 

The  OAREA  set  Identification  assigned  to  this  static  load  vector     (Integer  >  0) 

Load   set  identification  number  of  a   set  of  static   load  cards    (Any  card   that  may  be 
referenced  by  the  L3A0  Case  Control    card)      (Integer  >  0  or  blanic) 

Set   Identification   of   a    thermal    load    set    (Any  card   that  may   be   referenced   by    the 
TEMP(L0AO)  Case  Control    card)      (Integer  >  0  or  blanic) 


Remarks:      1.       The    above   cards   will    not  be  used   unless    selected   in   the   Case   Control    Deck   with    a 
L0AOSET  card. 

2.  This   card    is    available  only    in    suoerelements   and    1n   Solution    Sequences    25,    27,    30 
and  31. 

3.  The    number   of    static    load    vectors    created    for   each    suoerelement    is    the   number   of 
unique  OAREA   IDs  on  all    LSEQ  cards   in  the  bulk  data. 

4.  The  OAREA   identification   assigned   to   the   static   load   vectors  may   be   r«»fe-°nc3d   by 
RL3A01,   RL0AO2,   TL3A01   and  TL0AO2  cards. 

5.  Element  data  recovery  for  thermal    loads  is  not  currently   implemented   in  dynamics. 

6.  OAREA    set   identification   numbers    should  be   unique   with   respect   to    all    static    load 
set  identification  numbers. 
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Input  Data  Card   TABLED1   Oynamic  Load  Tabular  Function,  Form  I 

Description:    Defines  a  tabular  function  for  use  In  generating  frequency-dependent  and  time- 
dependent  dynamic  loads 

Format  and  Example: 


1 

2 

3 

4 

5 

5 

7 

9 

9 

10 

TA8LED1 

id    \^><Z1Z><Z^><CI^><C^ 

TABLED  1 

32 

1                                         '                     i 

[♦ABC 

xl        |        *l 

x2 

n 

x3 

*3 

*4           n 

*8C 

-3.0 

5.9 

2.0 

5.6 

3.0 

5.5 

ENOT 

etc- 


F1eld 


ID 


*1-*1 


Contents 
Table   Identification   number    (Integer  >  0) 
Tabular  entries    (Real) 


Remarks:      1.       The  x.,   must   be   in  either  ascending  or  descending  order  but    not   both. 

2.  Jumps   between   two  points   (x,    ■  xf+.)   are  allowed,    but   not   at   the  end   points. 

3.  At    least  two  entries  must   be  present. 

4.  Any  x-y  entry  may  be  Ignored  by  placing  the  3CD  string  "SKI?"  in  either  of  the  two 
fields  used   for  that  entry. 

5.  The  end  of  the  table  Is  Indicated  by  the  existence  of  the  3CD  string  "ENOT"  in 
either  of  the  two  fields  following  the  last  entry.  An  error  is  detected  if  any 
continuation   cards    follow   the   card  containing  the   end-of-table   flag    "ENOT". 

6.  Each  TA8LED1  mnemonic  Infers  the  use  of  a  specific  algorithm.  For  TA8LED1  type 
tables,   this   algorithm  1s 


•yt(x) 


7. 


where  X  1s  Input  to  the  table  ajid  i  is  returned.  The  table  look-up  yT(x),  x  -  (, 
1s  performed  using  linear  interpolation  within  the  table  and  linear  extrapolation 
outside  the  table  using  the  last  two  end  points  at  the  appropriate  table  end.  At 
jump  points  the  average  y-r(x)  Is  used.  There  are  no  error  returns  'rom  this  table 
look-up  procedure.  „ 

Linear  extrapolation  is  not  used  for  Fourier  Transform  methods.  The  function  is 
zero  outside  the  range. 
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Input  Data   CaM        [STEPNL        "  ~w    Step   tnd   Nonlinear    Sc    .'    .-    Contra    | 

0— Crf Pt low:      Provides    (MraaetHc    control*    Md    cata    for   ion    'n?Jr   Tmnimnt    fell    .< -.    | 


'ormat    and   Exano'e: 


1 

2 

3 

4 

5 

5 

7 

3 

? 

10 

T3TEPNL 

ID 

fOT 

9 

■<3 

HETH0O  | 

(STEP 

VUlTEfl 

33'./ 

rSTEPNL 

250 

^ 

: 

■jTfl 

; 

£ 

-  •. 

£?:;j 

£?:? 

E?3« 

MMOI  / 

^;:i 

FSTR£$2 

m.i 

.31 

a 

5 

.32 

Pi  »i  ^  3or».  ent t 

10  Identification   -j-oe*-   [Integer  >  0) 

10T  Numcer   of    time   steps    of    <a'.ie   3"'         '-*.?-:»'    /    3 

0T  Time   increment    ''ea;    /   3.0) 

•<0  3<1p    'actor   'or   output    'eve-y   "ifltn    step    -- '      :e    :a'»0    'or    v.:;.-  "»yt-    ,    '. 

"E~J0O  Metr<od    for   controlling    Iterat'ons    a"d  tancent   -at--/    .'.-.-.'.■■■.'.      333    «    -."3    jr    "3~"'- 

(No  default ) 

<3TEP  luraer  pf   steps   oe'ore  a  -at—*   .peat?    [Integer  ^_  3       3e'ij   t    *   3 

MAX:"'?  Limit   on  total    Iterations    'or   eac~   v-w    ' -c- <**-<'- 1      i-teoe"   •    3       '.r'-.    *.    i    3 

C3NV  r'aos     to    se:ect     ponve'-ze-ce    tests     '333     »    J,     5,     «,     or    My     cone inat    .'     p'     9ie 

I  etters    or  olanic ) 

3;3J  ,£?S?  ,£?sw     Er^or  tolerance  'iraits   'or  determining   convergence     '»*     >  '..'. 
[Defaults  ■  1.3E-2,   l.OE-3,   L.OE-4,   respective 

*AXD17  _imit    on    tre    -"jree-    ~f    diverging    iterations    et    eacJi    tine    and    m         er^  --    *.•-*» 

steps   [Integer  >  3)    [Oefaalt  *  I 

*axqn  Maximum   mjreer    of    3uas l-Hewt an    correction    lectori    -.p    :•»    :•:"»':    on    tha    data     -    - 

(Integer  >_  3)    [Default   «   23 

*axis  -'Jx-aum     nurce-     of      Line      3ea-p-       terations      for      .'.'.'•      notion        Integer   -   '. 

[Oefaalt  »  2)  " 

FSTBCSS  r~act'on    of    e"ecv  m    st'^ss   fnj    jsec    :p    pet?'-'*    the    Hibincrenent    i   EC     n    die 

naterial   rtwtlnes    '3.   <  Real    c  ..;      pefai  t  ■  '...'-.   .     3*e  Kenarfc    .3. 
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NASTRAN   DATA   OEC1C 


TSTEPNL    (Cont.) 


Remarks:  1.  The  TSTEPNL  Bulk  Data  card  Is  selected  by  the  Case  Control  card  TSTEPNL  -  ID.  Each 
suOcase  (residual  suoerelement  solutions  only)  requires  a  TSTEPNL  card  and  either 
aopHed  loads  via  TL0AD1  data  or  Initial  values  from  a  previous  subcase.  Multiple 
subcases  are  assumed  to  occur  sequentially  1n  time  with  the  initial  values  of  time 
and  displacement  conditions  of  each  subcase  defined  by  the  end  conditions  of  Che 
previous  subcase. 

2.   OT,  SOT,  HQ   art   the  time  increments,  At,  the  total  numOer  of  steps  Nt,  and  the  out- 
put interval,  nn,  respectively.  The  value  of  time,  Tn,  at  each  Increment  n  will  be 


tn  -  tQ  ♦  nat     n  »  0,  1 N, 


where  tg  1s  the  Initial  value  of  time  which  corresponds  to  the  last  steo  from  the 
previous  subcase. 

3.  For  printing  and  plotting  the  solution,  data  will  be  output  at  steps  n  »  0,  n^, 
2ng,....N.  The  Case  Control  card  0TIME  may  also  be  jsed  to  control  the  output 
points. 

4.  The  basic  solution  method  Is  the  "Newmarx  9eta"  transient  numerical  method  along 
with  the  "Newton-Raphson"  method  for  controlling  nonlinear  effects.  Equilibrium 
solutions  with  respect  to  lnertial,  damping,  elastic,  and  nonlinear  loads  ir<* 
obtained  at  each  time  steo  much  the  same  as  In  the  statics  case. 

The  KMETH0O  options  provide  the  user  with  three  ootional  nonlinear  iteration 
methods.  These  ire: 

AUTJ  provides   for  automatic   stiffness  matrix  updates   to   improve  poor 
convergence.  The  KSTEP  value  Is  Ignored. 

TSTEP  provides  for  a  matrix  update  at  every  KSTEPtn  increment  of  time 
(typi  cal ly  a  large  number). 

5.  The  stiffness  matrix  1s  always  updated  for  a  new  subcase  or  restart.  Only  one 
update  1s  allowed  for  each  time  step  and  the  Iteration  process  returns  :o  the 
previous  time  steo  after  in   update. 

6.  MAXITER  defines  the  number  of  static  nonlinear  load  iterations  (corrections)  to  be 
performed  at  each  time  step.  At  least  two  iterations  ire  necessary  to  determine 
the  convergence  of  the  displacements  (see  below). 

7.  The  test  flags  (U,P,'«0  and  the  error  limits  (EPSU,  etc.)  define  the  convergence 
criteria  (U  ■  Displacement  error  test,  P  »  Load  equilibrium  error  test,  W  *  'Wonc 
error  test).  If  the  internally  calculated  error  fractions  ir»  less  than  the 
requested  limits,  the  Iteration  stoos,  the  results  are  processed,  and  the  program 
continues  to  the  next  time  step  or  subcase. 
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TSTEPNl  (Cont.) 


8.  If  io  tests  are  requested  or  no  convergence  occurs,  the  numDer  of  iterations  is 
limited  to  MAXITER.  The  results  will  be  processed  and  the  program  «ill  continue 
with  the  next  time  step. 

9.  MAXQN  provides  the  control  over  diverging  solutions.  4t  each  time  step,  the 
solution  1s  allowed  to  diverge  until,  at  the  MAXDIV  iteration,  a  matrix  uodate  ;s 

requested.   If  the  average  convergence  parameter  T  Indicates  divergence  for  *^AXOIV 

successive  time  steps,  the  transient  solution  stoos  and  exits  for  data  recovery. 

10.  The  number  of  subincrements  1n  the  material  routines  (elastoplastlc  and  creep)  in 
determined  such,  that  the  subincrement  size  1s  approximately  FSTRESS  *  7. 

FSTRESS  Is  also  used  to  establish  a  tolerance  for  elastoplastlc  material  in  the 
material  routine,  i.e.. 


6   in  yield  function  <  FSTRESS  •  yield  stress 


If  the  limit  Is  exceeded  at  the  converging  state,  the  program  «i  1 1  EXIT  *i  th  a 
fatal  error  message.  Otherwise  the  stress  state  Is  adjusted  to  the  current  yield 
surface,  resulting  In  5  »  Q. 
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I  APPENDIX  C 


OUTPUT  OF  MSC/NASTRAN 
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